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ABSTRACT 


This  report  describes  the  results  of  experimental  and  theoretical 
investigations  of  the  mechanisms  by  which  the  core  flow  combustion  process  in 
coaxial,  single  inlet,  dump  type  ramjet  engines  drives  longitudinal 
combustion  instabilities.  To  this  end,  the  behavior  of  V-shaped  flames, 
similar  to  those  often  occurring  in  ramjet  combustors,  stabilized  in 
longitudinal  acoustic  fields  has  been  studied.  The  presence  of  burning 
vortical  structures  is  observed  in  the  flame  region.  These  structures  appear 
at  frequencies  close  to  the  first  natural  acoustic  frequency  of  the  combustor 
and  are  believed  to  be  connected  with  a  shear  layer  type  of  instability  of 
the  flame.  Experiments  conducted  show  that  the  unsteady  combustion  in  these 
structures  is  capable  of  driving  the  acoustics  at  the  fundamental  acoustic 
mode  frequency.  With  increase  in  fuel  air  ratio,  a  spontaneous  instability 
involving  the  fundamental  mode  is  observed  and  explained  in  terms  of 
increased  driving  associated  with  the  higher,  unsteady  heat  release  rates. 
The  results  of  experiments  conducted  with  external  acoustic  excitation  of  the 
flame  at  different  frequencies  are  also  reported  and  confirm  the  idea  that 
the  vortical  structures  arise  due  to  a  fluid  mechanical  instability  of  the 
flame.  It  is  shown  that  the  interactions  between  the  vortical  structures  and 
the  system's  acoustic  field  affect  the  heat  release  rates  from  the  flame  and 
provide  a  mechanism  for  the  driving  of  longitudinal  mode  instabilities. 


Chapter  I 


INTRODUCTION 


The  primary  objective  of  the  research  program  described  in  this  report 


was  the  determination  of  the  characteristics  of  the  driving  by  the  combustion 


process  which  is  responsible  for  the  occurrence  of  low  frequency  longitudinal 


instabilities  in  liquid  fueled,  single  inlet,  coaxial  dump  type  ramjet 


engines  (see  Fig.  1).  Recent  developments  of  such  ramjet  engines,  as  part  of 


integral  rocket-ramjet  propulsion  systems1  ,  have  been  hindered  by  the 


occurrences  of  destructive  combustion  instabilities 


These  instabilities 


are  characterized  by  either  low  frequency  (i.e.,  rumble)  or  high  frequency 
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(i.e.,  screech)  pressure  and  velocity  oscillations  ’  .  The  low  frequency 


rumbTe  is  in  the  range  of  50-500  Hertz  and  it  is  generally  characterized  by 


longitudinal  acoustic  oscillations  in  the  inlet  section  and  the 
combustor*1  ’  ^ 14  .,4-^,.  n 


or  a  combination  of  a  longitudinal  acoustic  oscillation  in 


the  inlet  section  and  a  bulk-mode  type  oscillation  in  the  combustor  .  In 


contrast,  the  high  frequency  screech  occurs  when  one  of  the  tangential 


acoustic  modes  of  the  combustor  is  excited  (e.g.,  see  Table  I  in  Ref. 5).  Of 


the  two,  the  low  frequency,  longitudinal  instability  presents  a  more  serious 


problem  because  the  interaction  of  the  resulting  pressure  oscillations  with 


the  inlet  shock  system  may  result  in  inlet  unstarting  and  loss  of  engine 


performance.  In  addition,  this  type  of  instability  may  result  in  excessive 


vibrational  loads  on  the  system.  On  the  other  hand  screech  type  instabilities 


result  in  an  increase  in  heat  transfer  rates  to  engine  components  which  may 


shorten  the  engine's  life  time  and  compromise  its  performance.  This  type  of 


instability  can,  however,  be  controlled  by  use  of  acoustic  liners,  unless  the 


throws* 


IDSTEeAM  AMD  FLUSH 
'ALL  fUZL  IOj£4TOet> 


engine  is  large  and  the  transverse  modes  are  characterized  by  low 
5 

frequencies  . 


While  both  types  of  instabilities  have  the  potential  for  severely 
hindering  the  development  of  present  and  future  ramjet  engine  programs,  of 
the  two,  the  low  frequency,  longitudinal  instability  is  potentially  more 
troublesome  because  of  the  lack  of  effective  means  for  its  suppression.  At 
present  there  is  no  proven  model  capable  of  predicting  the  conditions  under 
which  such  an  instability  may  occur  and  its  elimination  often  involves  a 
costly  and  time  consuming  trial  and  error  approach.  To  develop  a  more 
rational  approach  for  the  treatment  of  low  frequency  ramjet  instabilities, 
understanding  of  the  processes  responsible  for  their  occurrence  needed  to  be 
extended.  This  problem  was  addressed  in  the  research  program  described  in 
this  report. 


To  date,  combustion  instabilities  have  appeared  in  almost  every  type  of 

propulsion  system,  including  liquid16  and  solid16  rocket  motors,  air 
17  18 

breathing  engines  ,  ramjets  and  so  on.  In  the  majority  of  cases  the 

combustion  process  provides  the  energy  reauired  for  the  excitation  and 

maintenance  of  the  observed  oscillations.  These  oscillations  also  experience 

energy  losses  due  to  viscous  dissipation,  heat  conduction,  acoustic  energy 
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radiation  and  convection  through  the  exhaust  nozzle  ’  .  During  an 
instability,  the  amplitude  of  the  combustor  oscillation  grows  in  time  as  long 
as  the  energy  added  to  the  oscillation  per  cycle  is  larger  than  the  energy 
lost  per  cycle  due  to  the  above  mentioned  loss  processes.  As  the  amplitude 
grows,  some  or  all  of  the  energy  gain  and  loss  processes  become  amplitude 
dependent  and  an  amplitude  is  reached  at  which  the  energy  added  to  the 


oscillation  per  cycle  equals  the  energy  lost  per  cycle.  When  a  balance 
between  the  energy  addition  and  removal  processes  is  established,  the 
amplitude  of  the  oscillation  remains  constant  as  long  as  the  above  mentioned 
balance  is  not  disturbed. 

An  unstable  engine  can  be  stabilized  by  either  decreasing  the  energy 
supplied  to  the  oscillation  by  the  combustion  process  (or  any  other  process) 
or  increasing  the  disturbance  losses  or  by  doing  both.  To  accomplish  either 
of  these  steps,  a  thorough  understanding  of  the  processes  in  question  is 
required. 
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According  to  Rayleigh's  criterion  ’  ,  driving  of  the  combustor's 
oscillations  by  the  combustion  process  will  occur  if 

T 

E  = j  j p  Q  dtdV  >  0  (1) 

V  o 

In  the  above  integral,  E,p,Q,t,T,  and  V  represent  the  total  energy  added  to 
the  system  per  cycle,  the  pressure  oscillation,  the  oscillatory  heat  addition 
by  the  combustion  process,  time,  period  of  the  oscillation  and  volume  of  the 
system,  respectively.  The  integration  is  performed  over  the  volume  of 
interest  to  account  for  all  possible  energy  addition  processes.  When  p  and  Q 
are  in  phase,  the  above  integral  is  positive  and  vice  versa.  Also,  p  and  Q 
might  be  in  phase  in  some  portion  of  V  and  out  of  phase  in  the  remainder  of 
V.  In  such  a  case,  the  sign  of  E  will  depend  upon  the  relative  contributions 
of  the  "driving"  and  "damping" ( i .e. ,  sections  where  p  and  Q  are  out  of  phase) 
fractions  of  V.  Combustion  instability  will  occur  when  the  integral  E, 


Eq.(l),  is  larger  than  some  quantity  L  which  describes  the  total  losses 
experienced  by  the  system's  oscillation  per  cycle. 

In  coaxial,  single  inlet,  dump  type  ramjet  combustors  (see  Fig.  1), 

liquid  fuel  is  injected  into  the  air  flow  in  the  inlet  section  upstream  of 

the  dump  plane.  The  liquid  fuel  drops  evaporate  and  mix  with  the  air  as  they 

both  move  towards  the  combustor.  The  flow  expands  as  it  enters  the  combustor 

and  a  recirculation  region,  which  serves  as  a  flame  holder,  is  established 

near  the  step.  While  it  is  not  clear  what  fraction  of  the  fuel  burns  in  the 

recirculation  zone,  it  is  expected  that  most  of  the  fuel  burns  in  the  flame 

around  the  central,  jet-like,  core  flow  (  see  Fig.  1)  which  disappears  within 

a  few  inlet  section  diameters  downstream  of  the  dump  plane.  The  resulting 

ramjet  flow  field  is  extremely  complex  involving  turbulent  mixing,  a 

recirculating  flow  region,  a  jet  flow  region  and  finite  width  combustion 

zones.  In  addition,  shear  layers  are  formed  by  the  flow  at  the  dump  as  well 

as  by  the  flame  itself  which  separates  the  relatively  low  speed  unburned  flow 

from  the  relatively  higher  speed  burned  (reacted)  flow.  These  shear  layers, 

via  a  fluid  mechanical  instability  mechanism,  may  roll  up  into  vortices  as 
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has  been  confirmed  in  the  present  ’  ramjet  studies  described  in  this 
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report  as  well  as  by  other  research  groups  ’  .  Hence,  it  is  not  surprising 

that,  to  date,  attempts  to  model  steady  state  ramjet  flows  have  only  met  with 
limited  success  and  they  have  involved  some  empiricism  as  well  as  the  latest 
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advances  in  computational  modelling  techniques  .  This  complexity 
increases  considerably  when  one  attempts  to  model  the  unsteady  operation  of  a 
ramjet  engine,  as  is  required  in  combustion  instability  analyses. 


The  question  of  whether  periodic  heat  release  by  the  combustion  process 
in  a  ramjet  occurs  and  if  so  whether  it  occurs  in  phase  with  the  pressure 
oscillation  as  required  by  Rayleigh's  criterion  for  combustion  instability 
(i.e.,  see  Eq.(l))  depends  upon  the  response  of  the  combustion  process  to 
local  pressure  and  velocity  oscillations.  For  example,  local  pressure  and 
velocity  oscillations  may  initiate  a  fluid  dynamical  instability  of  the  flame 
causing  it  to  roll  up  into  vortices.  If  this  vortex  shedding  becomes 

periodic,  the  unsteady  combustion  in  the  vortices  would  also  become  periodic. 
If  the  resulting  periodic  heat  release  rate  is  in  phase  with  the  local 

pressure  oscillation,  an  amplification  of  the  oscillation  will  occur  as 
indicated  by  Eq.  (1)  above.  Consequently,  it  is  important  to  develop  an 
understanding  of  the  combustion  process  response  to  local  oscillations  (or 

disturbances)  as  such  an  understanding  may  lead  to  means  for  reducing  or 

minimizing  combustion  instabilities  in  ramjets. 


The  current  research  program  has  been  concerned  with  this  understanding 
of  the  basic  mechanisms  of  ramjet  instability.  The  evidence  to  date  supports 
the  idea  that  a  major  culprit  causing  ramjet  instabilities  is  the  unsteady 
heat  release  from  large  scale  coherent  vortical  structures  that  exist  in  the 
core  flame  region  of  the  ramjet  combustor.  A  sequence  of  frames  taken  by 
means  of  high  speed  shadowgraphy  in  the  experimental  set  up  developed  at 
Georgia  Tech  to  study  ramjet  instabilities  is  shown  in  Fig. 2.  The  flame  is 
stabilized  in  a  rectangular  combustor  by  a  heated, thin  nichrome  wire  (see 
Fig. 3).  The  vortical  structures  are  clearly  visible.  Note  also  that  the 
vortical  structures  arise  symmetrically  on  the  top  and  lower  flame  fronts  and 
cause  finite  width  combustion  zones.  As  will  be  discussed  in  the  next 
chapter,  it  is  believed  that  these  vortical  structures  arise  by  means  of  a 


shear  layer  type  of  instability  mechanism  of  the  flame  region.  These  vortices 
are  convected  along  the  flame  front  at  approximately  the  local  flow  velocity. 
The  phase  of  the  unsteady  heat  release  from  the  combustion  in  these  vortices 
is  related  to  their  convection  velocity.  As  the  flow  in  the  combustor  section 
is,  generally,  of  low  Mach  number  the  convection  velocity  is  small  compared 
to  the  local  speed  of  sound  which  determines  the  phase  of  the  acoustic 
oscillations.  Hence,  the  phase  of  the  acoustic  pressure  oscillation  changes 
slowly  as  compared  to  the  phase  of  the  unsteady  heat  release.  Thus,  over 
certain  regions  of  the  flame  the  unsteady  heat  release  and  the  pressure 
oscillation  are  in  phase  while  over  other  regions  they  are  out  of  phase. 
Details  are  provided  in  the  next  chapter  and  in  Appendices  B  and  C.  The 
important  fact  is  that  when  any  significant  pressure  oscillations  (100  dB  or 
more)  are  measured  in  the  set  up  it  is  found  that  the  "driving"  or  the  in 
phase  part  of  the  flame  exceeds  the  "damping"  or  out-of-phase  part  of  the 
flame.  Therefore,  in  accordance  with  Eq.  (1),  the  unsteady  heat  release  from 
the  combustion  in  the  vortical  structures  inputs  energy  into  the  acoustic 
oscillations.  However,  this  by  itself  does  not  mean  that  the  sole  cause  of 
the  observed  pressure  oscillations  is  the  unsteady  heat  release  from  the 
vortical  structures. 

To  determine  whether  the  observed  pressure  fluctuations  are  indeed  caused 
by  the  unsteady  heat  release  from  the  vortical  structures  a  theory  has  also 
been  developed  for  predicting  the  pressure  oscillations  from  the  unsteady 
heat  release  from  the  flame.  Briefly,  the  pressure  spectrum  is  constructed 
from  the  heat  release  rate  by  means  of  the  following  equation  (derived  in 
Appendix  D) 


Here,  Spp  is  the  pressure  spectrum,  Spp  is  the  unsteady  heat  release  spectrum 
and  G  is  the  relevant  Green's  function  which  describes  how  an  unsteady  heat 
source  gets  "converted"  into  pressure.  Now,  if  the  predicted  spectrum  from 
Eq.  (2)  agrees  with  the  experimentally  measured  pressure  spectrum  in  the  set 
up  then  it  may  be  concluded  that  indeed  the  unsteady  heat  release  from  the 
vortical  structures  is  the  main  cause  of  the  observed  acoustic  oscillations. 
That  this  is  the  case  is  confirmed  by  comparisons  between  the  two  spectra  (an 
example  is  presented  in  Fig.  10). 

Further  details  of  the  investigations  carried  out  under  the  research 
program  are  described  in  this  report.  In  Chapter  II,  the  major  results  of  the 
experimental  and  theoretical  efforts  are  documented.  The  short  comings  of 
recent  theoretical  models  of  ramjet  instability  are  also  considered  there. 
Conclusions  and  suggestions  for  further  research  are  presented  in  Chapter 
III.  The  report  closes  with  a  number  of  appendices. 


REVIEW  OF  RAMJET  INSTABILITY  STUDIES  AT  GEORGIA  TECH 


The  primary  objectives  of  the  ramjet  combustion  instability  research 

program  at  Georgia  Tech  described  herein  were:  (i)  to  determine  the 

characteristics  of  the  combustion  process  in  ramjet  combustors  responsible 

for  the  occurrence  of  low  frequency,  longitudinal  instabilities  in  liquid 

fueled,  coaxial,  dump  type  ramjet  engines  and  (ii)  to  investigate  the 

validity  of  previously  developed  linear  combustion  instability  models,  such 
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as  the  model  developed  by  Yang  and  Culick  . 

To  meet  these  objectives  the  set  up  shown  in  Fig.  3  was  developed  (see 

also  Appendices  B  and  C).  It  consists  of  an  inlet,  combustor  and  exhaust 

sections.  The  inlet  contains  a  movable  injector  through  which  a  mixture  of 

propane  and  air  is  introduced  into  the  combustor.  The  flame  is  stabilized  in 
2 

the  7.5  x  5  cm  combustor  section  on  a  0.8  mm  diameter  nichrome  wire.  The 
wire  is  attached  to  the  combustor  windows  at  half  combustor  height  and  it  is 
heated  electrically  to  improve  the  steadiness  of  the  stabilized  v-shaped 
flame.  The  exhaust  section  is  equipped  with  two  acoustic  drivers  which  are 
used  to  excite  a  standing  acoustic  wave  of  desired  amplitude  and  frequency  in 
the  set  up.  This  allows  the  study  of  the  effects  of  known,  external 
excitation  on  the  flame.  In  addition,  the  movable  injector  provides  a 
capability  for  "placing"  the  stabilizing  wire  (and  hence,  the  flame)  on  any 
part  of  the  excited  standing  wave;  that  is,  at  a  pressure  maximum,  minimum  or 
in  between.  The  maximum  length,  1,  available  between  the  injector  face  and 
the  exhaust  section  is  3  meters.  The  flow  approaching  the  stabilizing  wire  is 


parallel  and  uniform  and  disturbances  are  damped  out  by  means  of  a  fine  wire 

mesh  grid  located  8  cm.  upstream  of  the  stabilizing  wire.  This  wire  mesh  also 

serves  as  a  flame  arrestor  in  case  of  flash  back.  The  combustor  walls  are 

water  cooled  enabling  wall  mounted  pressure  transducers  to  be  used  to  monitor 

the  acoustic  pressure  field.  Capabilities  for  time  and  space  resolved 

measurements  of  CH  species  radiation  from  the  flame  have  also  been  developed. 
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The  concentrations  of  these  species  are  a  measure  of  the  reaction  rate  and, 
thus,  the  heat  release  rate  and  are  useful  in  describing  the  unsteady 
combustion  field. 

The  acquired  pressure  and  radiation  signals  are  digitized  and  stored  in 
computer  memory  prior  to  their  fourier  analyses.  The  developed  fourier 
analysis  program  yields  both  auto  spectra  and  cross  spectra  of  the  signals. 
This  facilitates  determination  of  relationships  between  the  pressure  and  the 
unsteady  heat  release  including  their  relative  phases.  Also,  capabilities  for 
identifying  specific  frequency  components  and  their  amplitude  and  phase 
behavior  have  been  developed. 

The  acoustic  nature  of  the  set  up  was  first  investigated  under  cold  flow 

conditions.  It  was  determined  that  the  injector,  whose  face  is  made  of 

sintered  stainless  steel,  could  be  well  approximated  as  an  acoustically  rigid 

surface.  Thus,  since  the  exhaust  end  of  the  set  up  is  open,  the  developed  set 

up  behaves  as  an  acoustically  closed-open  organ  pipe.  The  natural  acoustic 
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frequencies  of  this  system  are,  therefore,  given  by  the  quarter  wave  mode  , 
the  three  quarter  wave  mode  and  so  on. 


The  behavior  of  the  flame  has  been  investigated  under  conditions  of  (i) 
no  driving  by  the  acoustic  drivers,  (ii)  with  driving  provided  by  the 
acoustic  drivers  at  different  frequencies  and  (iii)  in  a  shortened  version  of 
the  setup  which  was  considered  due  to  reasons  set  forth  shortly.  The  behavior 
of  the  flame  under  each  of  these  conditions  is  described  below. 

(i)  Results  Obtained  with  No  External  Driving:  Spontaneous  instabilities 
involving  the  set  up's  fundamental  acoustic  mode  (i.e.,  the  quarter  wave 
mode)  have  been  observed  in  certain  fuel/air  ratio  ranges.  When  the  mixture 
composition  is  close  to  the  lean  stability  limit,  the  flame  is  thin  and  can 
be  readily  stabilized  on  the  heated  wire.  However,  as  the  fuel  fraction  is 
increased,  instability  sets  in  as  indicated  by  the  wall  mounted  pressure 
transducers  which  record  an  increase  in  amplitude  of  the  fundamental  acoustic 
mode.  With  further  increase  in  the  fuel  fraction  (while  the  mixture  remains 
in  the  fuel  lean  range)  the  flame  flashes  back  from  the  wire  to  the  screen 
and  the  sound  pressure  levels  in  the  combustor  approach  140-150  db.  As  the 
mixture  becomes  fuel  rich,  not  all  the  fuel  can  be  burned  within  the 
combustor  and  part  of  the  fuel  burns  outside  the  set  up.  When  the  mixture 
approaches  the  rich  flammability  limit  the  flame  is  blown  back  to  the  wire 
where  it  stabilizes  again.  Under  these  conditions  the  instability  subsides 
and  the  flame  is  again  stable. 

As  observations  of  the  flame  indicated  that  an  instability  (as  recorded 
by  the  increasing  pressure  levels  with  increase  in  fuel  air  ratio)  set  in 
while  the  flame  was  still  held  at  the  wire,  experiments  were  undertaken  under 
these  conditions  to  investigate  the  cause  of  the  instability.  These  were 
investigated  by  use  of  CH  radiation  measurements.  By  a  suitable  arrangement 
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of  the  optical  elements,  the  radiation  from  vertical  2  mm  wide  strips  of  the 
flame  were  measured.  The  optics  could  be  moved  axially  so  that  the  radiation 
from  different  parts  of  the  flame  could  be  measured.  Figure  4  presents  the 
axial  dependences  of  the  radiation  amplitudes,  along  the  flame  front  (the 
wire  corresponds  to  x=0)  for  three  different  frequencies  of  40  Hz,  80  Hz  and 
160  Hz.  Note  the  sharp  increase  in  amplitude  of  all  three  frequencies  just 
downstream  of  the  wire.  The  first  natural  acoustic  frequency  is  80  Hz  and  it 
is  clear  that  the  signal  at  this  frequency  dominates  the  local  unsteady 
radiation  (and,  thus,  also  the  unsteady  heat  release). 

To  link  this  behavior  to  a  shear  layer  type  of  instability,  some  of  the 

features  of  shear  layer  instabilities  will  now  be  considered.  Shear  layer 
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instabilities  have  primarily  been  investigated  for  non  reacting  flows  ,  in 
particular  for  incompressible  flows  and  this  discussion  will  .therefore,  be 
limited  to  such  flows.  However,  as  such  instabilities  are  primarily  fluid 
mechanical  in  nature,  it  is  expected  that  their  major  features  will  also  be 
present  even  under  reacting  flow  conditions.  When  a  shear  layer  instability 
sets  in,  the  amplitudes  of  the  disturbance  (e.g.,  the  disturbance  velocity) 
increases  exponentially.  Also,  the  instability  is  confined  to  a  limited 
frequency  range  (i.e.,  the  Strouhal  number  range).  Disturbances  at 
frequencies  in  this  limited  range  are  amplified  whereas  disturbances  at  other 
frequencies  are  not  and  may,  instead,  be  damped.  A  typical  plot  of  the  growth 
parameter,  which  is  related  to  the  disturbance  growth  rate,  as  a  function  of 
the  relevant  Sti ouhal  number  is  shown  in  Fig.  5  taken  from  Ref.  34.  Note  that 
the  growth  parameter  is  non  zero  only  in  a  limited  frequency  range.  The 
frequency  at  which  the  growth  parameter  maximizes  is  the  most  unstable 
frequency. 
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Figure  5.  Typical  Variations  of  the  Nondimens ional 
Growth  Rate  with  Strouhal  Number  (taken 
from  Ref.  34). 


Looking  at  the  radiation  measurements  from  this  point  of  view,  the  sharp 
increase  in  the  radiation  measurements  just  downstream  of  the  wire  could 
indicate  a  similar  instability  behavior  of  the  flame  shear  layer.  More 
evidence  to  this  effect  is  obtained  by  noting  the  relative  amplitudes  of  the 
different  frequency  components.  While  the  signal  at  the  first  natural 
frequency  of  80  Hz  dominates,  the  signal  at  40  Hz  is  also  significant. 
However,  the  signals  at  frequencies  higher  than  80  Hz  start  diminishing  in 
magnitude.  For  example,  the  signal  at  160  Hz  is  much  smaller  than  both  the  80 
and  40  Hz  components.  It  was  found  that  signals  at  higher  frequencies 
(greater  than  160  Hz  )  rapidly  diminished  in  magnitude.  Thus,  it  is  possible 
that  the  flame  in  the  set  up  exhibits  a  shear  layer  type  of  instability  which 
is  confined  to  frequencies  up  to  approximately,  160  Hz. 

Additional  supporting  evidence  is  obtained  by  considering  the  phase  of 
the  radiation  signal  with  respect  to  the  pressure  oscillations.  Figure  6 
presents  the  phase  of  the  radiation  signal  at  the  first  natural  frequency  of 
80  Hz  as  a  function  of  the  axial  location  along  the  flame.  Note  the  smooth 
monotonic  variation  of  the  phase.  The  smoothness  of  the  phase  curve  indicates 
that  the  disturbances  at  80  Hz  are  well  correlated  and  coherent  over  the 
entire  flame  as  would  be  expected  if  they  were  the  result  of  a  shear  layer 
type  instability.  It  was  found  that  such  smooth  phase  curves  could  only  be 
obtained  for  the  frequencies  under  approximately  160  Hz  which,  as  estimated 
earlier,  was  the  extent  of  the  frequency  range  for  the  shear  layer  type  of 
instability  of  the  flame.  Although  it  will  be  discussed  under  the  results 
obtained  with  external  driving,  it  may  be  pointed  out  that  disturbances 
imposed  upon  the  flame  by  the  acoustic  drivers  in  the  frequency  range  of 


iation  Signal  (at  80  Hz)  with  Respect  to  Pressure 


160-500  Hz  were  rapidly  damped  out  by  the  flame.  This  result  also  confirms 
the  limited  frequency  range  of  the  flame  shear  layer  instability. 

The  phase  variation  of  the  low  frequency  components  may  be  related  to 
the  phase  velocity  of  the  disturbances.  This  is  done  by  noting  that  a  phase 
change  of  2rr  radians  (i.e.,  360°)  occurs  over  a  distance  that  may  be  regarded 
as  one  disturbance  wavelength.  The  convection  velocity  uc  is  then  given  by 
the  product  of  the  relevant  frequency  with  the  corresponding  disturbance 
wavelength,  provided  they  remain  constant.  In  the  more  general  case,  when 
these  quantities  vary  over  the  flame,  it  is  more  appropriate  to  consider  the 
convection  velocity  between  neighboring  points  along  the  flame.  In  such  a 
case  the  convection  velocity  (u  )..  between  two  points  x.  and  x.  on  the  flame 

w  l  J  J 

is  given  by 

2nf[x  -xj 

(u  )  *  ~ "  ^ 

ciJ  V9i  (3) 

where  9.  and  e.  are  the  phases  at  x.  and  x.,  respectively. 
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The  convection  velocity  at  80  Hz,  derived  from  the  above  equation  is 
plotted  in  Fig.  7  as  a  function  of  the  axial  distance  along  the  flame.  The 
value  at  the  flame  stabilizing  wire  is  about  1.5  m/sec  which  is  equal  to  the 
unburned  mixture  flow  velocity.  Its  value  at  the  downstream  end  of  the  flame 
is  about  6  m/sec  which  is  equal  to  the  burned  flow  velocity  calculated  on  the 
basis  of  the  density  drop  across  the  flame.  Thus,  the  disturbance  at  80  Hz  is 
convected  along  the  flame  at  approximately  the  average  of  the  local  unburned 
and  burned  ( i  .e.  .reacted)  flow  velocities.  This  fact  also  implies  that  the 
disturbances  in  the  flame  are  fluid  mechanical  in  nature.  If  they  had  been 


purely  acoustic,  their  convection  or  propagation  velocity  would  have  been  the 
local  speed  of  sound  and  not  the  local  flow  speed.  The  other  disturbances  in 
the  frequency  range  up  to  160  Hz  also  behave  similarly,  thus  demonstrating 
their  fluid  mechanical  nature. 

These  fluid  mechanical  disturbances  have  been  identified  as  burning 
vortical  structures  that  are  shed  from  the  flame  holding  wire  by  means  of 
high  speed  (6000  frames/sec)  shadowgraphy  and  schlieren  photography  performed 
on  the  flame.  It  was  found  that  the  frequencies  of  shedding  favored 
predominantly  frequencies  around  the  first  natural  frequency  of  the  set  up. 
This  was  found  to  be  true  even  when  external  driving  by  means  of  the  acoustic 
drivers  was  provided  as  will  be  described  shortly  (see  also  Fig.  2). 

Thus,  the  results  of  these  experiments  demonstrate  that  the  flame 
region  is  fluid  dynamically  sensitive  in  a  limited  frequency  range 
(approximately  up  to  160  Hz  in  the  developed  set  up)  and  is  involved  in  a 
shear  layer  type  of  instability  in  this  frequency  range.  This  results  in  the 
formation  of  vortical  disturbances  which  are  convected  along  the  flame  front. 
The  unsteady  combustion  in  these  structures  favors  frequencies  in  the  same 
frequency  range  that  the  flame  exhibits  the  shear  layer  type  of  instability. 
The  question  now  arises  whether  the  unsteady  combustion  is  responsible  for 
the  pressure  oscillations  that  lead  to  the  observed  spontaneous  instability. 

To  answer  this  question  it  was  checked  whether  the  unsteady  heat  release 
in  the  sensitive  frequency  range  satisfied  Rayleigh's  criterion  for  flame 
driving  of  the  acoustics  (i.e.,  see  Eq.  (1)).  For  convenience,  Eq.  (1)  may  be 
rewritten  in  the  frequency  domain  in  the  following  form 
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S  =  the  cross  spectrum  between  the  unsteady  heat  release  and  the 
PQ 

pressure  oscillations  at  a  given  frequency  and 

=  the  corresponding  phase  between  the  pressure  and  the  heat 

release 

and  Vx,  is  the  volume  of  the  flame  region. 
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The  integrand  of  the  above  equation  is  plotted  in  Fig.  8  for  the  first 
natural  frequency  of  80  Hz.  As  the  flame  length  is  much  smaller  than  the 
acoustic  wavelength  at  80  Hz  the  pressure  is  constant  over  the  flame  region. 
Thus,  the  oscillatory  nature  of  the  integrand  is  related  to  the  wavelength  of 
the  disturbance  or  vortical  structures  at  80  Hz.  Since  the  integral  was  found 
to  be  positive,  these  results  indicate  that  the  unsteady  heat  release  at  80 
Hz  drives  the  acoustic  oscillations  at  that  frequency.  The  Rayleigh  criterion 
was  also  found  to  be  satisfied  at  other  low  frequency  components  indicating 
that  the  flame  also  inputs  energy  into  the  acoustic  oscillations  (i.e., 
drives)  at  these  frequencies. 

Satisfying  the  Rayleigh  criterion  does  not  exclude,  however,  the 
presence  of  other  important  acoustic  energy  sources  in  the  experimental  set 
up.  For  example,  impingement  of  the  convected  vortices  against  a  solid 
surface  (e.g.,  the  duct  walls)  could  also  contribute  to  the  excitation  of 


duct  acoustic  oscillations  as  has  been  suggested  to  occur  in  segmented  solid 
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propellant  rocket  motors  .  Thus,  it  becomes  necessary  to  predict  the 
acoustic  pressure  spectrum  by  considering  the  unsteady  heat  release  as  the 
primary  acoustic  source  and  comparing  this  prediction  with  measured  pressure 
spectra.  If  satisfactory  agreement  is  obtained,  then  the  unsteady  heat 
release  may  be  taken  to  be  the  primary  acoustic  source. 

A  theory  was,  therefore,  developed  to  predict  the  pressure  spectrum  from 
the  unsteady  heat  release.  According  to  this  theory,  valid  for  low 
frequencies  for  which  the  flame  zone  may  be  taken  to  be  compact  compared  to 
the  acoustic  wavelength  (see  derivation  in  Appendix  D). 

SPP  •  fl  1°!2  S<,q  (2) 

P 

Here,  Spp  is  the  acoustic  pressure  spectrum,  Sqp  is  the  net  or  integrated 
heat  release  spectrum  from  the  entire  flame  and  G  is  the  relevant  Green's 
function  which  describes  how  an  unsteady  heat  source  excites  acoustic 
pressure  oscillations. 

The  integrated  heat  release  spectrum  was  obtained  by  measuring  the  CH 
radiation  from  the  entire  flame.  The  result  is  shown  in  Fig.  9.  Note  that  now 
the  radiation  at  the  first  natural  frequency  does  not  dominate  the  spectrum 
as  it  did  for  the  local  heat  release  at  different  locations  along  the  flame. 
The  reason  for  this  is  that  the  flame  length  is  equivalent  to  several 
vortical  wavelengths  as  is  evident  from  Fig.  8.  This  means  that  when  the 
signal  is  integrated  over  the  entire  flame,  a  lot  of  cancellation  occurs  from 
the  positive  and  negative  portions  of  the  unsteady  heat  release.  When  lower 


frequencies  are  considered,  due  to  their  greater  wavelength,  this 
cancellation  is  smaller.  As  a  result  ,  it  is  the  frequencies  lower  than  80  Hz 
and,  in  particular,  frequencies  under  50  Hz  that  dominate  the  spectrum. 

A  typical  comparison  of  the  predicted  and  measured  pressure  spectra  is 
shown  in  Fig.  10.  The  comparison  is  excellent  and,  thus,  confirms  that  the 
unsteady  heat  release  is  the  primary  acoustic  source  in  the  set  up.  It  may  be 
noted  that  although  there  are  substantial  levels  of  pressure  oscillations  at 
frequencies  lower  than  80  Hz,  the  dominant  peak  in  the  pressure  spectrum 
occurs  at  the  first  natural  frequency  of  80  Hz.  This  takes  place  because  the 
Green's  function  acts  like  a  filter  and  favors  frequencies  close  to  the 
natural  acoustic  frequencies  of  the  duct  (see  Appendix  D).  Thus,  although  the 
source  strengths  at  the  lower  frequencies  are  higher  (see  Fig.  10)  the 
pressure  oscillations  favor  the  natural  acoustic  frequencies  of  the 
experimental  setup. 

It  was  also  found  that  with  increase  in  the  fuel/air  ratio  the  magnitude 
of  the  unsteady  heat  release  (the  source  for  the  pressure  oscillations)  also 
increases.  For  example,  Fig.  11  plots  the  measured  radiation  signal  as  a 
function  of  the  burned  gas  temperature  for  the  natural  frequency  of  80  Hz. 
The  measurement  was  made  3  cm  downstream  of  the  flame  holding  wire  and  the 
temperatures  are  normalized  by  the  temperature  at  which  flash  back  occurs. 
Note  the  sharp  increase  in  the  signal  as  the  burned  gas  temperature  (fuel 
fraction)  increases. 

Hence,  the  spontaneous  instability  observed  in  the  set  up  can  be 
explained.  The  shear  layer  formed  by  the  flame  undergoes  a  fluid  mechanical 
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type  of  instability  causing  the  shedding  of  vortices  from  the  flame  holding 
wire.  The  shear  layer  type  of  instability  is  confined  to  a  certain  frequency 
range.  The  resulting  unsteady  combustion  drives  the  acoustic  oscillations  in 
the  set  up.  With  increase  in  the  fuel  fraction  the  unsteady  heat  release 
rates  increase.  The  consequent  increase  in  the  flame  driving  then  leads  to 
the  observed  instability. 

(ii)  Results  Obtained  with  External  Driving:  Experiments  were  conducted  to 
study  the  effects  of  external  excitation,  provided  by  the  acoustic  drivers, 
on  the  flame  behavior.  These  reinforced  the  concept  of  a  shear  layer  type  of 
instability  of  the  flame  which  is  confined  to  a  limited  frequency  range.  The 
driving  was  provided  at  different  frequencies  in  the  range  of  160-500  Hz. 

It  was  found  in  all  these  cases  that  the  flame  could  be  divided  into  two 
regions,  that  is,  (a)  a  region  close  to  the  flame  holding  wire,  called  the 
near  wake  region,  where  the  dominant  radiation  signal  is  at  the  driving 
frequency  and  (b)  a  region  downstream  of  the  near  wake  region,  called  the  far 
wake  region  where  the  lower  frequency  components,  in  particular  the  signal  at 
the  first  natural  frequency,  dominate.  Thus,  even  under  external  driving,  the 
flame  region  exerts  a  dominating  influence  in  causing  the  decay  of  the  non 
sensitive  or  stable  frequency  disturbances.  An  example  is  shown  in  Fig.  12 
where  the  near  wake  and  far  wake  regions  for  driving  at  385  Hz  are  shown.  The 
near  wake  and  far  wake  behavior  was  found  to  be  independent  of  the  location 
of  the  flame  holding  wire  on  the  excited  standing  wave;  i.e.,  at  a  pressure 
maximum,  a  pressure  minimum  or  in  between. 


385  Hz 


These  results  were  confirmed  by  high  speed  schlieren  and  shadow 
photography  (6000  frames/sec)  of  the  flame  under  external  excitation.  A 
sequence  of  frames  taken  from  the  shadow  graph  movie  is  shown  in  Fig.  2.  The 
vortical  structures  were  still  found  to  formed  around  the  first  natural 
frequency  of  the  set  up.  Visually,  the  only  effect  of  the  external  driving 
was  an  oscillation  of  the  vortical  structures,  as  they  were  convected  along 
the  flame  front,  at  the  driving  frequency. 


Results  Obtained  in  a  Shorter  Version  of  the  Setup: 


As  the 


experiments  described  above  indicated  that  the  unstable  Strouhal  number  range 
for  the  flow  conditions  utilized  in  the  setup  and  the  first  natural  frequency 
of  the  setup  overlapped,  it  was  felt  that  a  situation  where  this  did  not 
occur  should  also  be  investigated.  Hence,  a  shorter  version  of  the  setup  was 
considered.  The  exhaust  section,  downstream  of  the  combustor  section,  was 
removed.  This  shortening  resulted  in  an  increase  in  the  first  natural 
frequency  of  the  setup.  Depending  upon  the  injector  location  upstream  of  the 
flame  holding  wire,  the  first  natural  frequency  of  the  setup  could  be 
increased  to  about  200  Hz  which  was  outside  the  sensitive  Strouhal  number 
range. 


As  before,  radiation  measurements  at  different  locations  along  the  flame 
as  well  as  int'  rated  radiation  measurements  from  the  entire  flame  were 
obtained.  A  strong  radiation  component  at  the  first  natural  frequency  was 
found  to  be  present  only  very  close  to  the  flame  holding  wire.  Downstream  of 
this  region,  the  signal  at  the  first  natural  frequency  decayed  and  the  low 
frequency  components  dominated  (see  Fig.  13).  In  this  respect,  the  first 
natural  frequency  behaved  like  an  externally  driven  frequency  outside  the 
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sensitive  Strouhal  number  range,  in  the  longer  setup  as  described  in  (ii). 
However,  it  was  found  that  unlike  in  the  longer  setup,  the  low  frequency 
components  in  the  sensitive  Strouhal  number  range  were  not  well  correlated 
over  the  flame  front.  In  particular,  a  well  defined  convection  velocity  could 
not  be  obtained  at  these  low  frequencies.  For  the  signal  at  the  first  natural 
frequency  also,  a  well  defined  convection  velocity  did  not  exist  over  the 
entire  flame.  A  convection  velocity  could  be  calculated  at  the  first  natural 
frequency,  f  ,  in  this  case  only  in  the  near  portion.  This  is  shown  in  Fig. 
14,  for  f  =  180  Hz.  However,  this  velocity  seems  to  bear  little,  if  any, 
relation  to  the  mean  flow  velocity  as  in  the  longer  setup.  In  the  far  wake 
region,  moreover,  the  convection  velocity  is  not  calculable  to  any  reasonable 
precision  due  to  uncorrelated  changes  in  the  phase  behavior  of  the  signal. 
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A  possible  explanation  for  this  behavior  is  that  some  amount  of  acoustic 
feedback  is  necessary  to  correlate  disturbances  over  the  entire  flame  even 
for  those  frequencies  which  lie  in  the  Strouhal  number  range  for  the  shear 
layer  instability.  As  shown  in  Appendix  0,  the  amount  of  acoustic  feedback 
from  the  flame  at  a  given  frequency  is  dependent  upon 

a)  the  source  strength;  i.e.,  the  integrated  unsteady  heat  release  from  the 
entire  flame  at  the  given  frequency  and 
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b)  how  far  is  this  given  frequency  from  the  nearest  natural  acoustic 
frequency. 


CONVECTION  VELOCITY 
NOT  CALCULABLE 
(SHORTENED  SET  UP) 


Natural  Frequencies  in  the  Long 


Thus,  when  the  first  natural  acoustic  frequency  becomes  removed  from  the 
shear  layer  instability  regime,  the  acoustic  feedback  at  these  frequencies 
will  decrease  according  to  b)  above. 


To  test  the  effect  of  acoustic  feedback,  the  shortened  setup  was  excited 
by  the  acoustic  drivers  at  the  natural  frequency  of  f  =  180  Hz.  Figures  15 
and  16  show  the  effect  of  this  driving.  Note  the  sharp  increase  in  the 
amplitude  of  the  radiation  at  180  Hz  in  the  near  wake  (Fig.  15).  It  is  also 
seen  that  the  extent  of  the  near  wake  has  increased  compared  to  the  no 
driving  case  (compare  Figures  13  and  15).  Also,  from  Fig.  16,  it  is  clear 
that  with  the  driving,  a  well  defined  convection  velocity  which  is  close  to 
the  velocity  of  the  mean  flow,  exists  over  the  entire  flame.  This  confirms 
the  importance  of  an  adequate  amount  of  acoustic  feedback  to  correlate 
unsteady  disturbances  over  the  flame. 


li cations  for  Ramjet  Combustion  Instability  Models: 


The  results  of  the  experiments  described  above  shed  light  on  the 
deficiencies  of  current  ramjet  combustion  instability  models.  These  aspects 
will  now  be  discussed. 


Referring  to  Fig.l  it  is  clear  that  the  ramjet  flow  field  is  extremely 
complex  involving  turbulent  mixing,  a  recirculating  flow  region  in  the 
vicinity  of  the  dump  plane,  a  jet  flow  region,  shear  layers  and  finite  width 
combustion  zones.  Thus,  even  the  modeling  of  only  steady  state  ramjet 
combustor  flow  fields  presents  the  investigator  with  a  myriad  of  challenging 
problems  which  are  related,  in  one  way  or  another,  to  the  lack  of  complete 


V  *.’  * .*  V  ■/  ■/  *,*  V  »  r  *  "  « ■ 


35 


180  Hz(f  ) 


Axial  Distance  from  Flame 
(Shortened  Set  Up). 


understanding  of  the  complex  ramjet  combustion  and  flow  processes.  These 
difficulties  increase  when  attempts  are  made  to  model  the  unsteady  operation 
of  a  ramjet  engine  as  is  required  in  combustion  instability  analyses. 

Notwithstanding  these  complexities,  efforts  have  been  made  by  some 
31  36  37 

investigators  ’  ’  to  develop  ramjet  combustion  instability  models.  Here, 
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the  model  developed  by  Yang  and  Culick  will  be  discussed  to  assess  the 
applicability  of  such  models.  This  model  appears  to  reasonably  predict  the 
acoustic  wave  structure  in  the  combustor  but  fails  to  predict  stability 
limits.  In  the  Yang  and  Culick  model  (to  be  denoted  henceforth  by  YCM) 
considerable  simplification  is  achieved  by  assuming  that  the  amplitude  of  the 
disturbance  is  small,  which  reduces  the  mathematics  to  the  solution  of  a 
system  of  linear  differential  equations.  As  a  result  of  these  simplifications 
the  model  can  only  hope  to  predict  the  operating  conditions  under  which 
spontaneous  combustion  instability  may  occur  and  the  wave  structure  of  the 
small  amplitude  oscillations.  The  model  requires  that  a  solution  (or  a 
description)  of  the  steady  ramjet  combustor  flow  field  be  available  because 
the  steady  state  solutions  appear  in  the  coefficients  of  the  equation  which 
describe  the  behavior  of  the  instability. 

To  make  the  problem  tractable,  a  number  of  assumptions  are  made  in  the 
YCM.  Notable  among  these  is  the  assumption  that  all  of  the  combustion  in  the 
ramjet  combustor  takes  place  inside  a  thin  flame  sheet.  The  model  assumes 
that  when  an  oscillation  occurs  in  the  combustor,  the  flame  sheet  gets 
distorted  by  the  perturbations,  but  remains  infinitesimally  thin.  The  flame 
area  is  a  measure  of  the  amount  of  combustibles  consumed  by  the  flame  and 
periodic  variation  in  the  flame  area  may  result  in  a  periodic  heat  release 


rate  which,  in  turn,  may  drive  acoustic  pressure  oscillations  when  the  proper 
phase  relationship  between  the  heat  release  rate  and  the  local  pressure 
oscillations  is  satisfied.  This  is  the  primary  acoustic  driving  mechanism 
considered  in  the  YCM. 

While  the  experiments  conducted  during  this  investigation  do  not  negate 
the  idea  that  the  oscillatory  heat  release  from  the  flame  arises  due  to 
periodic  variations  in  flame  area  they  do  demonstrate  that  any  variations  in 
the  flame  surface  are  caused  by  the  presence  of  vortical  structures  in  the 
combustion  zone.  There  is  no  consideration  of  vortical  structures  in  the 
flame  zone  in  the  YCM.  Also,  as  is  clear  from  the  shadowgraph  film  sequence 
in  Fig  2  the  combustion  zone  is  not  a  zero  thickness  flame  sheet  so  that  the 
concept  of  flame  surface  area  must  be  carefully  defined. 

The  longitudinal  acoustic  instability  of  the  combustor  is  linked  to  the 
unsteady  combustion  rate  at  an  axial  acoustic  natural  frequency  in  the 
vortical  structures.  The  unsteady  combustion  rate  favors  those  frequencies  to 
which  the  flame  shear  layer  is  fluid  mechanically  unstable.  Thus,  if  an  axial 
acoustic  natural  frequency  of  the  combustor  overlaps  the  unstable  frequency 
range  of  the  flame  shear  layer  high  rates  of  unsteady  combustion  may  be 
obtained  at  this  frequency  which,  in  turn,  could  lead  to  high  levels  of 
acoustic  oscillations  or  the  onset  of  longitudinal  combustion  instability. 

Thus,  it  is  clear  that  the  acoustic  or  longitudinal  combustion 
instability  of  the  ramjet  must  be  considered  in  conjunction  with  the  fluid 
mechanical  type  of  instability  of  the  flame  shear  layer.  So  far,  fluid 
dynamicists  have  focussed  on  shear  layer  instabilities  in  non  reacting  flows 


with  much  of  the  work  being  done  in  incompressible  flows.  A  fluid  dynamic 
instability  theory  for  a  chemically  reacting  shear  layer  does  not  exist  at 
present.  Such  a  theory  is  needed  to  predict  the  frequency  range  in  which  the 
burning  vortical  structures  appear  and  to  understand  their  response  to 
changes  in  the  flow  and  chemical  composition  of  the  combustible  mixture. 


Chapter  III 


CONCLUSIONS 


From  the  results  of  the  experimental  studies  described  earlier,  the 
importance  of  vortex  shedding  and  the  unsteady  combustion  in  these  vortices 
in  causing  and  maintaining  ramjet  combustor  instabilities  is  evident.  The 
unsteady  combustion  in  the  shed  vortices  feeds  energy  into  the  combustor 
acoustics  if  the  following  inequality  (i.e.,  Rayleigh's  criterion)  is 
satisfied  (Eq.  (4)) 
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where  S  =  cross  spectrum  between  pressure  and  the  unsteady  heat  release 
PQ 

6pq  =  phase  between  the  pressure  and  the  unsteady  heat  release  and 

V,,  =  volume  of  the  flame  region 

t  i  ame 

As  |  Spq  |  >  0,  the  sign  of  the  above  integral  is  largely  determined  by 
the  variation  of  the  phase  angle  8pq  over  the  flame  region.  Thus,  over  some 
regions,  the  flame  may  be  driving  while  over  other  regions  the  flame  may  be 
damping. 


The  spectrum  of  the  pressure  fluctuations  Spp  is  related  to  the  spectrum 
of  the  unsteady  heat  release  S  ,  at  least  for  compact  flames,  by  the 
following  equation  (Eq.  (2)) 


where  G  is  the  relevant  Green's  function  and  describes  how  an  unsteady  heat 
source  excites  pressure  oscillations. 

Since  G  is  maximized  at  the  natural  acoustic  frequencies  of  the 
combustor  (Appendix  D),  then  Eq.  (2)  indicates  that  if  high  values  of 
unsteady  heat  release  occurs  at  a  natural  acoustic  frequency  of  the 
combustor,  large  pressure  fluctuations  will  also  occur  at  that  frequency.  If 
vortex  shedding  occurs  at  a  natural  Trequency  of  the  combustor,  then,  large 
values  of  Sqq  may  be  obtained  at  this  frequency. 

From  these  considerations  it  is  clear  that  to  suppress  combustion 
instabilities  one  or  more  of  the  following  approaches  should  be  undertaken: 

(i)  vortex  shedding  (  at  a  natural  frequency)  must  be  suppressed  and/or 

(ii)  the  phase  angle  6pq  must  be  so  controlled  that  the  inequality,  Eq.  (2) 
is  not  satisfied  and/or 

(iii)  the  magnitude  of  Sqq  must  be  reduced. 

Further  research  is  needed  to  understand  how  best  these  approaches  may 
be  incorporated  in  dump  type  ramjet  combustors.  For  example,  the  shear  layer 
instability  of  the  flame  should  be  investigated  to  determine  the  conditions 
under  which  vortex  shedding  may  occur  at  an  acoustic  natural  frequency  of  the 
combustor.  This  could  lead  to  guidelines  for  designing  the  flame  holder 


region  to  suppress  vortex  shedding  at  an  acoustic  natural  frequency.  The 
fundamental  aspects  of  the  interaction  between  the  geometry  of  the  flame  zone 
and  its  fluid  mechanics  should  also  be  studied  to  gain  understanding  of  how 
the  magnitude  and  phase  of  the  unsteady  heat  release  may  be  controlled. 

In  conclusion,  it  is  believed  that  this  research  program  has  been 
successful  in  identifying  one,  if  not  the  major,  source  of  longitudinal 
instabilities  in  ramjet  like  combustors;  that  is,  the  modification  of 
combustion  rates  in  the  ramjet  flame  zone  by  vortical  structures  arising  as  a 
result  of  the  fluid  mechanical  instability  of  the  flame. 
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ABSTRACT 

Coaxial,  dump  type  ramjet  combustors  are  often 
prone  to  combustion  instability  problems  that  can 
seriously  affect  their  performance.  This  study  is 
concerned  with  the  highly  detrimental  low  frequen¬ 
cy  instabilities  which  represent  one  of  the  most 
common  types  of  the  encountered  instabilities.  In 
particular,  the  coupling  between  the  cone-like 
flames  which  are  stabilized  at  the  entrance  of  the 
dump  combustor  and  longitudinal  acoustic  fields  is 
studied  in  an  experimental  set  up  specifically 
developed  for  this  purpose.  It  consists  of  a 
rectangular  tube  having  an  injector  plate  at  one 
end  and  acoustic  drivers  at  the  opposite  end.  A 
combustible  mixture  enters  through  the  injector 
and  a  conical  flame  is  stabilized  on  a  thin  wire 
attached  to  the  side  walls  of  the  tube.  The 
behavior  of  this  flame  in  various  longitudinal 
acoustic  fields  generated  by  the  acoustic  drivers 
Is  studied.  A  variety  of  experimental  techniques, 
including  high  speed  Schlieren  and  shadow  pnotog- 
raphy,  CH  radiation  emission  and  acoustic  pressure 
measurements  are  used.  Investigations  were 
undertaken  at  two  experimental  set  up  lengths.  A 
strong  coupling  between  the  flame  response  and 
certain  natural  acoustic  frequencies  is  observed 
in  the  longer  set  up.  For  the  shorter  set  up,  the 
coupling  is  seen  to  be  intense  when  the  flame 
stabilization  region  is  located  near  an  acoustic 
pressure  minimum.  These  observations  shed  new 
light  on  flame/acoustic  mode  interactions  in  a 
ramject  combustor-like  environment  and  the  limita¬ 
tions  of  current  state  of  the  art  theoretical 
models  of  this  phenomenon. 


INTRODUCTION 

Recent  developments  of  coaxial  dump  type 
ramjet  engines  have  been  hindered  by  occurrences 
of  destructive  combusion  instabilities  '  .  These 

instabilities  are  characterized  by  either  low 
frequency  (i.e.  rumble)  or  high  frequency  (i.e. 
screech)  pressure  and  velocity  oscillations.  The 
low  frequency  rumble  is  in  the  range  of  100-500 
Hertz  and  is  generally  characterized  by 
longitudinal  acoustic  oscillations  in  the  inlet 
section  and  the  combustor.  In  contrast,  the  high 
frequency  screech  occurs  when  one  of  the 
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tangential  acoustic  modes  of  the  combustor  is 
excited.  In  the  present  study  the  low  frequency 
type  of  instability  is  considered.  This  type  of 
instability  can  interfere  with  the  shock  system  in 
the  ramjet  inlet  resulting  in  loss  of  engine 
performance  as  well  as  imposing  excessive 
vibrational  loads  on  the  system. 

Combustion  instabilities  have  appeared  in  a 
variety  of  propulsion  systems  including  liquid  and 
solid  propellant  rocket  motors,  air  breathing 
engines,  ramjets  and  so  on.  In  the  majority  of 
cases  the  combustion  process  provides  the  energy 
required  for  the  excitation  and  maintenance  of  the 
observed  oscillations.  According  to  Rayleigh's 
criterion  '  ,  driving  of  the  pressure  oscilla¬ 

tions  by  the  combustion  process  occurs  if  the 
oscillating  heat  release  from  the  combustion 
process  is  in  phase  with  the  pressure  oscilla¬ 
tions.  Under  these  conditions,  energy  is  added  to 
the  oscillations  during  each  cycle.  However,  the 
acoustic  field  also  loses  energy  due  to  viscous 
dissipation,  heat  conduction  and  the  convection 
and  radiation  of  acoustic  energy  through  the 
exhaust  nozzle.  During  an  instability,  tne 
amplitude  of  the  combustor  oscillation  grows  in 
time  as  long  as  the  energy  added  to  the  oscilla¬ 
tions  per  cycle  is  larger  than  the  energy  lost  per 
cycle  due  to  the  above  mentioned  loss  processes. 

As  the  amplitude  grows,  some  or  all  of  the  energy 
gain  and  loss  processes  become  amplitude  dependent 
and  an  amplitude  is  reached  at  which  the  energy 
added  to  the  oscillation  per  cycle  equals  the 
energy  lost  per  cycle.  When  this  condition  is 
established,  the  amplitude  of  the  oscillations 
remains  constant. 


To  date,  there  has  been  a  paucity  of 
experimental  data  relating  to  the  coupling  between 
the  combustion  processes  and  acoustic  oscillations 
in  a  ramjet  engine  like  environment.  Recent 
studies  include  those  of  Davis  and  Heitor  et 
al  .  With  the  help  of  high  soeed  photography 
Davis  categorized  two  modes  of  low  frequency 
instability  in  dump  type  combustors.  In  tne  first 
mode,  the  entire  combustion  zone  underwent  a 
cyclic  oscillation  while  in  the  second  mode 
regular  shedding  of  hot  spots  from  the 
recirculation  zone  at  the  dump  plane  was  observed. 
Heitor,  et  al  measured  the  frequency  and  strength 
of  combustion  induced  oscillations  for  premixed 
methane  air  flames  stabilized  on  baffles  located 
on  the  axis  of  a  pipe.  They  found  that 
instability  involving  the  first  natural  mode 
(quarter  wave  type)  occurred  for  a  wide  range  of 
fuel  to  air  ratios.  More  recently.  Smith  and 
Zukoski  '  '  studied  combustion  Instability 
Involving  a  flame  stabilized  behind  a  rearward 
facing  step.  They  report  that  under  conditions  of 
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Instability  larga  vortical  structures  Mere  formed 
downstream  of  the  step  at  an  acoustic  resonant 
frequency.  These  were  subsequently  convected  away 
from  the  step  and  the  resulting  unsteady 
combustion  was  found  to  feed  energy  1ntogthe 
acoustic  oscillations.  Keller,  et  al 
Identified  three  modes  of  Instability  In  a  similar 
set  up  in  an  effort  to  understand  the  mechanisms 
responsible  for  flashback.  They  attributed  this 
phenomenon  to  the  action  of  vortices  In  the 
recirculation  zone  of  the  step.  These  studies 
Indicate  that  the  fluid  mechanics  of  the  flame 
stabilization  region  plays  an  important  role  in 
combustion  Instability. 

Theoretical  efforts  to  model  combustion 
instabilities  In  dump  type  combustors  may  be 
broadly  classified  into  two  types.  The  first  type 
assumes  the  flame  zone  to  be  compact  relative  to 
the  acoustic  wavelength.  The  governing  flow 
equations  on  either  side  of  the  flame  zone  are 
matched  by  making  several  assumptions  about  the 
combustion  process.  The  analyses, gf  Abouseif  et 
al  and  Dowling  and  Bloxsldge  •  '  exemplify 
this  approach.  The  second  type. of  analysis,  as 
developed  by  Yang  and  Cullck  '  ,  Investigates  the 

combustion  region  in  more  detail.  They  used  an 
Idealized  description  of  the  flow  field  inside  a 
ramjet.  Their  model  assumes  the  combustion  to 
take  place  in  an  Infinitely  thin  flame  sheet  which 
Is  distorted  by  the  presence  of  pressure  and 
velocity  oscillations.  Using  integral  forms  of 
the  conservation  equations  Yang  and  Cullck  ob¬ 
tained  solutions  for  both  the  combustor  acoustic 
wave  structure  and  the  complex  frequency  of 
oscillation  whose  imaginary  part  determines 
whether  an  Instability  will  occur  or  not. 

Unfortunately,  neither  of  these  approaches 
have  been  totally  successful.  Whereas  proponents 
of  the  first  approach  claim  to  be  able  to  predict 
instability  no  understanding  of  the  coupling 
between  the  acoustic  field  and  the  unsteady 
combustion  process  Is  obtained  since  the  combus¬ 
tion  zone  acoustic  wave  structure  Is  ignored.  The 
Yang-Cull ck  model  on  the  other  hand  appears  to 
reasonably  predict  the  acoustic  pressure  structure 
In  the  combustor  but  fails  to  predict  stability 
limits.  Thus,  there  Is  a  need  to  understand  In 
more  detail  the  complex  processes  that  link  the 
acoustic  field  to  the  unsteady  combustion  field. 
The  present  study  focuses  on  this  interaction. 

EXPERIMENTAL  EFFORTS 

The  set  up  described  In  Fig.  1  has  been 
developed.  It  consists  of  an  Inlet,  combustor  and 
exhaust  sections.  The  inlet  contains  a  movable 
Injector  through  which  a  mixture  of  propane  and 
air  Is  Introduced  Into  the  combustor.  The  flame 
Is  stabilized  In  the  7.5  x  5  an  combustor  section 
an  a  0.8  mm  diameter  nlchrome  wire  chosen  to 
reduce  vortex  shedding  tendencies.  The  wire  Is 
attached  to  the  combustor  windows  at  half  combus¬ 
tor  height  and  It  Is  heated  electrically  to 
Improve  the  steadiness  of  the  stabilized  V-shaped 
flame.  The  exhaust  section  Is  equipped  with  two 
acoustic  drivers  which  are  used  to  excite  a 
standing  acoustic  wave  of  desired  amplitude  and 
frequency  In  the  set-up.  The  movable  Injector 
provides  a  capability  for  "placing"  the  stabiliz¬ 


ing  wire  on  any  part  of  an  excited  standing  wave; 
that  Is,  at  a  pressure  maximum,  minimum  or  in 
between.  The  maximum  length,  L  ,  available 
between  the  injector  face  and  tnexexhaust  section 
exit  plane  Is  3  meters.  The  flow  approaching  the 
stabilizing  wire  is  parallel  and  uniform  and 
disturbances  are  damped  out  by  means  of  a  fine 
wire  mesh  grid  located  8  cm.  upstream  of  the 
stabilizing  wire.  The  cold  flow  Reynolds  number 
is  kept  below  10,000  in  the  interest  of  maintain¬ 
ing  a  disturbance  free  flow.  The  combustor  walls 
are  water  cooled  enabling  wall  mounted  pressure 
transducers  to  be  used  to  monitor  the  acoustic 
pressure  field.  Capabilities  for  time  and  space 
resolved  measurements  of  CH  species  radiation  from 
the  flame  have  also  been  developed.  The  concen¬ 
trations  of  these. species  is  a  measure  of  the 
reaction  rate  '  ,  and,  hence,  the  heat  release 

rate  and  are  useful  in  describing  the  unsteady 
combustion  field. 

Experiments  have  been  performed  in  the  set  up 
described  above  and  in  a  modified  version  in  which 
the  exhaust  section  was  removed  and  the  acoustic 
driver  section  was  attached  directly  at  the 
combustor  exit.  The  consequent  shortening  of  the 
tube  results  in  higher  values  for  the  fundamental 
acoustic  frequencies  of  the  set  up.  Striking 
differences  were  observed  in  the  flame  response  to 
acoustic  excitation  in  the  two  set  ups.  These  are 
described  later. 

The  injector  plate  Is  made  of  sintered 
stainless  steel  and  It  acts,  approximately,  as  an 
acoustically  rigid  termination  so  that  the  acous¬ 
tic  behavior  of  the  experimental  set-up  approxi¬ 
mates  that  of  a  closed-open  organ  pipe.  During 
tests  In  both  set  ups,  spontaneous  instabilities 
involving  the  tube's  fundamental  acoustic  mode 
( 1 . e . ,  the  quarter  wave  mode)  have  been  observed 
in  certain  fuel/air  ratio  ranges.  When  the 
mixture  composition  Is  close  to  the  lean  stability 
limit,  the  flame  is  thin  and  can  be  readily 
stabilized  on  the  heated  wire.  However,  as  the 
fuel  fraction  in  increased,  instability  sets  in  as 
indicated  by  the  wall  mounted  pressure  transducers 
which  record  an  increase  in  amplitude  of  the 
fundamental  acoustic  mode.  With  further  increase 
in  the  fuel  fraction  (while  the  mixture  remains  in 
the  fuel  lean  range)  the  flame  flashes  back  from 
the  wire  to  the  screen  and  the  sound  pressure 
levels  in  the  combustor  approach  140-150  dB.  As 
the  mixture  becomes  fuel  rich,  not  all  the  fuel 
can  be  burned  within  the  combustor  and  part  of  the 
mixture  burns  outside  the  set-up.  When  the 
mixture  approaches  the  rich  flammability  limit  the 
flame  is  "blown  back"  to  the  wire  where  it  stabi¬ 
lizes  again.  Under  these  conditions  the  instabil¬ 
ity  subsides  and  the  flame  Is  again  stable  and 
thin. 

Observations  of  the  flame  Indicated  that  the 
Instability  originates  at  the  wire.  Thus,  experi¬ 
ments  conducted  to  date  have  been  carried  out  with 
the  flame  stabilized  at  the  wire.  These  included 
visualization  studies  with  high  speed  Schlieren 
and  shadow  photography,  CH  radiation  measurements 
from  the  flame  and  measurements  of  acoustic 
pressure  distributions  along  the  flame.  These 
have  been  carried  out  with  acoustic  excitation 
provided  by  the  drivers  at  different  frequencies. 


High  speed  (6000  frames/second)  Schlieren  and 
shadowgraph  films  were  taken  in  the  longer  tube 
with  acoustic  excitation  at  a  frequency  of  385  Hz 
and  acoustic  pressure  amplitudes  of  the  order  of 
130  dB.  By  varying  the  injector  location,  tests 
were  conducted  with  the  stabilizing  wire  placed  at 
a  pressure  maximum,  a  minimum  and  in  between.  The 
fuel  mass  fraction  was  2.2%  and  the  burned  gas 
temperature  as  measured  by  a  thermocouple  was  1325 
*K.  The  acoustic  pressure  along  the  flame  was 
simultaneously  monitored.  The  films  reveal  that 
there  is  a  periodic  shedding  of  flamelets  from  the 
wire.  As  one  flamelet  is  shed  another  is  formed 
at  its  tail  so  that  the  flame  is  continuous. 

These  flamelets  move  along  what  may  be  called, 
approximately,  the  steady  flame  front,  stretching 
as  they  proceed.  A  sequence  from  the  shadowgraph 
film  taken  with  the  wire  at  a  pressure  maximum  is 
presented  in  Fig.  2.  While  the  mixture  injection 
velocity  was  1.5m/sec,  the  convection  velocity  of 
the  flamelets  was  of  the  order  of  2.5m/sec. 
indicating  that  the  flamelets  are  convected  at  the 
local  flow  velocity.  The  surprising  fact  is  that, 
except  when  the  stabilizing  wire  was  at  a  pressure 
minimum  (i.e.  acoustic  velocity  maximum)  for  the 
driving  frequency  of  385  Hz.  the  shedding  was  at 
frequencies  corresponding  to  the  first  natural 
acoustic  mode  of  the  experimental  set-up  and  not 
at  the  driving  frequency.  However,  the  pressure 
amplitudes  of  the  first  natural  mode  were  at  least 
7db  below  those  of  the  acoustically  driven  signal 
at  385  Hz.  When  the  wire  was  moved  to  the  loca¬ 
tion  of  the  pressure  minimum  for  the  driven 
signal,  the  shedding  frequency  increased,  but 
still  remained  lower  than  the  driving  frequency. 
The  only  effect  characterized  by  the  driving 
frequency  in  all  the  investigated  cases  was  an 
observed  oscillation  of  the  flamelets,  as  they 
moved,  at  the  driving  frequency.  This  was  most 
evident  in  the  film  with  the  wire  located  at  the 
pressure  minimum  (i.e.,  a  velocity  antinode). 

To  investigate  this  phenomenon  further  CH 
radiation  measurements  were  carried  out  in  the 
long  tube.  By  a  suitable  arrangement  of  the 
optical  elements,  the  radiation  from  vertical  3  mm 
wide  strips  of  the  flame  were  made.  The  optics 
could  be  moved  axially  so  that  the  radiation  from 
different  parts  of  the  flame  could  be  measured. 
Initially,  the  measurements  were  carried  out  under 
conditions  similar  to  those  under  which  the  high 
speed  films  were  taken.  It  was  found  in  these 
cases  that  even  though  the  pressure  spectra  were 
dominated  by  the  contribution  of  the  385  Hz 
oscillation,  the  radiation  spectra  were  dominated 
by  the  frequencies  corresponding  to  the  first 
natural  frequency  for  each  case.  For  example,  in 
Fig.  3  the  amplitude  of  the  radiation  is  plotted 
as  a  function  of  the  axial  distance  from  the 
stabilizing  wire.  It  is  seen  that  the  magnitudes 
of  the  signals  at  the  first  natural  frequencies 
are  an  order  of  magnitude  larger  than  those  at  the 
driven  frequency.  Figure  A  describes  the 
frequency  spectrum  of  the  radiation  1  cm.  down¬ 
stream  of  the  wire  for  the  case  when  the  wire 
was  between  a  minimum  and  maximum  for  the  driving 
frequency.  It  is  seen  that  the  radiation  spectrum 
contains  peaks  at  frequencies  that  are  multiples 
of  the  first  natural  frequency.  For  the  ideal 
case  of  an  acoustically  closed-open  organ  pipe, 
which  the  system  approximates,  the  nth  longitudi¬ 


nal  acoustic  mode  is  f  =  (2n  +  1)  f  where  f  is 
the  frequency  of  the  f?rst  mode  (n  -0).  WitR  the 
addition  of  heat  due  to  combustion  and  the  result¬ 
ing  spatial  dependence  of  the  speed  of  sound  this 
relationship  becomes  only  approximate.  Therefore, 
it  should  be  noted  that  multiples  of  the  first 
natural  mode  frequency  which  are  seen  in  the 
radiation  spectrum  are  not,  in  general,  natural 
acoustic  frequencies  of  the  system.  The  pressure 
spectrum  for  this  case,  at  the  location  of  the 
wire  is  plotted  in  Fig.  5.  It  is  seen  to  be 
dominated  by  the  signal  at  385  Hz. 


The  phase  of  the  radiation  signal  at  the 
natural  frequencies  with  respect  to  the  pressure 
at  the  same  frequencies  for  the  above  cases  is 
plotted  in  Fig.  6.  A  mean  curve  has  been  drawn 
through  the  points  to  emphasize  the  nature  of  the 
data.  The  slope  of  the  phase  curve  is  inversely 
proportional  to  the  convection  velocity  of  the 
flamelets.  It  is  seen  to  decrease  with  increasing 
axial  distance  along  the  flame,  indicating  an 
increase  in  the  convection  velocity  with  distance. 
This  is  reasonable  because  of  the  increase  in 
temperature  that  occurs  along  the  flame  which 
results  in  an  acceleration  of  the  flow. 

Whereas  the  phase  of  the  radiation  signal 
changes  due  to  convection  by  the  local  flow,  the 
phase  of  the  acoustic  pressure  changes  according 
to  the  local  speed  of  sound  and  these  changes  are 
negligible  over  the  length  of  the  flame  for  the 
first  natural  frequency.  Thus,  over  certain 
regions  of  the  flame  the  radiation  and  the  acous¬ 
tic  pressure  are  in  phase  and  over  other  regions 
of  the  flame  they  are  out  of  phase.  It  is  inter¬ 
esting  to  note  that  the  maximum  of  the  amplitudes 
of  the  radiation  signal  occur,  approximately,  1-2 
cm  downstream  of  the  stabilizing  wire.  A  compari¬ 
son  of  Figures  3  and  6  shows  that  in  this  region 
there  is  a  component  of  the  radiation  at  the 
natural  frequency  which  is  in  phase  with  the 
pressure  thus  satisfying  Rayleigh's  condition  for 
driving. 

It  should  be  noted  that  the  above  tests  were 
carried  out  near  the  fuel  lean  flammability  limit 
with  the  fuel  mass  fraction  between  2.0  to  2.3%. 
When  the  fuel  fraction  is  increased  beyond  this 
level,  spontaneous  instability  occurs  as  described 
earlier.  This  spontaneous  instability  begins  even 
as  the  flame  is  attached  to  the  wire  and  with  a 
further  (small)  increase  in  the  fuel  mass  fraction 
the  flame  flashes  back  to  the  screen. 

Spectra  of  pressure  and  radiation  of  the 
unexcited  flame  are  shown  in  Figures  7  and  8  for 
different  fuel -air  ratios.  It  is  clearly  seen 
that  the  signal  strength  at  the  fundamental 
acoustic  frequency  (in  this  case  about  63  Hz) 
Increases  as  the  fuel-air  ratio  increases. 


The  above  discussed  data,  obtained  in  the 
“longer"  experimental  configurations.  Indicate 
that  there  is  a  strong  coupling  between  the 
radiation  oscillations,  which  are  a  measure  of  the 
unsteady  heat  release  rate  in  the  flame,  and  the 
fundamental  acoustic  natural  mode  frequency.  In 
an  effort  to  establish  whether  a  strong  coupling 


between  the  acoustics  and  the  flame  also  occurs  at 
other  frequencies,  tests  were  carried  out  with 
different  driven  frequencies  with  the  injector 
location  fixed.  The  first  four  acoustic  modes 
were  determined  to  be  at  63  Hz,  168  Hz,  282  Hz  and 
364  Hz,  respectively.  Since  the  drivers  are  not 
efficient  at  frequencies  below  200  Hz  where  a  lot 
of  the  acoustic  power  is  supplied  to  higher 
harmonics,  tests  were  run  at  frequencies  close  to 
the  third  and  fourth  natural  frequencies. 

When  the  experimental  set  up  was  driven  at 
the  fourth  mode  frequency  the  response  of  the 
flame  was  still  at  the  fundamental  frequency. 
However,  for  excitation  at  282  Hz  (the  third 
natural  frequency)  the  major  response  of  the  flame 
was  at  282  Hz.  However,  a  slight  shift  in  the 
driving  frequency  was  sufficient  to  change  this 
response.  For  example,  for  excitation  at  280  Hz, 
the  major  response  of  the  flame  was  no  longer 
around  the  third  natural  frequency.  This  is 
clearly  shown  in  Fig.  9  where  the  radiation 
signals  for  acoustic  excitations  at  282  Hz  and  280 
Hz  are  compared.  Thus,  a  response  at  the  same 
frequency  between  the  acoustics  of  the  system  and 
the  unsteady  behavior  of  the  flame  occurs,  for 
this  set  up,  only  at  a  few  selected  frequencies 
which  are  close  to  only  certain  natural  frequen¬ 
cies  of  the  system;  that  is,  in  the  above  case, 
the  first  and  third  and  not  at  all  the  natural 
frequencies. 

Experiments  have  also  been  conducted  with  the 
shorter  version  of  the  set  up  with  the  exhaust 
section  removed  and  the  acoustic  drivers  located 
directly  at  the  combustor  exit.  The  original 
reason  for  making  this  modification  was  to  get  the 
fundamental  natural  frequency  in  a  range  where  the 
acoustic  drivers  are  reasonably  efficient  thus 
enabling  the  study  of  the  flame  response  when 
excited  at  this  frequency.  Surprisingly,  it  was 
found  that  in  the  short  tube  the  flame  response 
did  net  always  occur  at  the  fundamental  natural 
frequency  when  the  excitation  was  at  different 
frequencies.  Instead  the  response  of  the  flame 
depended  strongly  upon  the  location  of  the  stabi¬ 
lizing  wire  on  the  generated  acoustic  wave  in  the 
range  200-500  Hz  in  which  the  tests  were  carried 
out.  For  example,  when  the  stabilizing  wire  was 
near  a  pressure  minimum  there  was  a  large  radia¬ 
tion  signal  at  the  driving  frequency.  On  the 
other  hand,  when  the  stabilizing  wire  was  near  a 
pressure  maximum  the  flame  response  was  not  at  the 
driving  frequency.  It  appears  that  the  acoustic 
velocity,  dominant  at  a  pressure  minimum  and  small 
at  a  pressure  maximum,  may  be  important  in  deter¬ 
mining  the  flame  response.  Experimental  work  on 
the  two  set  ups  is  continuing. 

Some  comments  on  stabilizing  the  flame  on  the 
heated  wire  are  also  in  order.  Originally,  the 
purpose  of  the  wire  was  to  both  ignite  and  stabi¬ 
lize  the  flame.  However,  igniting  the  mixture  by 
heating  the  wire  proved  to  be  quite  tricky. 

Hence,  ignition  was  achieved  by  use  of  a  propane 
torch  introduced  from  a  port  located  above  the 
wire.  Flames  stabilized  in  this  manner  on  the 
cold  wire  exhibited  a  great  degree  of  unsteadi¬ 
ness.  It  was  found  that  the  heating  of  the  wire 
greatly  reduced  this  tendency  and  also  helped  in 
avoiding  blow  off  of  the  flame.  It  did  not. 


however,  appear  to  have  any  effect  on  the  flame 
flashing  back  under  spontaneous  instability.  The 
wire  was  heated  by  an  alternating  current  at  60 
Hz.  The  fluctuating  component  of  the  wire  temper¬ 
ature,  due  to  ohmic  squaring,  was  therefore  at  120 
Hz.  No  signals  in  the  pressure  or  radiation 
spectra  were  found  at  this  frequency  that  could  be 
attributed  to  this  fluctuating  temperature  compo¬ 
nent  of  the  wire.  It  is  therefore  believed  that 
heating  of  the  wire  did  not  interfere  with  any  of 
the  phenomena  described  here. 

THEORETICAL  CONSIDERATIONS 

A, computer  program,  based  on  the  Yang-Culick 
model  ,  was  developed  to  solve  for  the  flow 
variables  in  the  combustion  zone.  Briefly,  it 
considers  the  flame  as  a  flow  discontinuity 
separating  the  unburned  and  burned  flows.  Conser¬ 
vation  equations  are  developed  and  solved  for  both 
the  upstream  and  downstream  regions  of  the  flame. 

A  kinematic  description  of  the  flame  front  and 
specification  of  the  flame  velocity  complete  the 
formulation.  The  assumption  of  small  perturba¬ 
tions  in  the  flow  variables  allows  the  split  of 
the  governing  equations  into  a  set  of  nonlinear 
differential  equations  for  the  steady,  time 
independent  part  and  a  set  of  linear  differential 
equations  for  the  unsteady  part.  Whereas  the 
Yang-Culick  model  obtains  the  natural  frequencies 
as  eigenvalues  of  the  time  dependent  equations, 
the  model  developed  for  the  present  study  assumes 
external  excitation  by  the  acoustic  drivers  and, 
hence,  a  known  frequency  of  oscillations. 

Observations  and  photography  of  the  flame 
indicate  that  under  steady  state  conditions,  the 
assumptions  of  the  Yang-Culick  model  appear  valid; 
that  is,  the  flame  is  thin  and  may  be  considered 
as  a  discontinuity.  Under  these  conditions  the 
flame  shape  can  be  reasonably  predicted.  An 
example  is  shown  in  Fig.  10.  The  experimental 
flame  shape  was  determined  by  direct  photography. 
The  theoretical  calculation  requires  the  values  of 
the  upstream  velocity  u  and  pressure,  the  flame 
velocity  W  (assumed  constant)  and  the  ratio  of  the 
temperature  across  the  flame.  The  flame  velocity 
W  was  measured  from  the  photograph  with  the 
knowledge  of  the  upstream  velocity  u  and  the 
flame  length.  In  all  cases,  it  was  round  that  the 
theory  underpredicted  the  flame  length.  This  may 
be  due  to  the  neglect  of  viscous  and  thermal 
conduction  effects  in  the  model.  These  effects 
become  important  near  the  wall  and  tend  to  lower 
the  flame  velocity  in  that  region  which,  in  turn, 
would  cause  lengthening  of  the  flame. 

The  unsteady  problem,  however,  poses  numerous 
hurdles.  The  Yang-Culick  model  assumes  that  in 
the  unsteady  case  the  flame  front  remains  thin 
though  it  undergoes  perturbations.  Viewing  of  the 
shadowgraph  films  and  the  phenomenon  of  the 
flamelet  shedding  makes  this  assumption  question¬ 
able.  Secondly,  the  distortion  of  the  thin  flame 
front  signifies  the  effect  of  the  combustion- 
acoustic  interaction  in  the  model.  Experimen¬ 
tally,  the  convection  of  the  flamelets  appears  to 
play  a  major  role  in  the  driving  of  the  tube's 
acoustic  modes.  The  assumption  of  a 
constant  flame  velocity  may  also  be  erroneous. 
Finally,  the  response  of  the  flame  at  a  frequency 
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other  than  the  driving  frequency  cannot  be  pre¬ 
dicted  by  a  linear  model. 

CONCLUSION 

The  response  of  wedge  shaped  flames,  similar 
to  those  often  occuring  in  ramjet  combustors,  to 
acoustic  excitation  has  been  studied  in  two  set 
ups  which  differed  by  the  length  of  their  exhaust 
sections.  For  the  longer  tube,  the  flame  responds 
strongly  to  imposed  excitation  only  at  certain 
preferred  frequencies.  These  frequencies  are 
close  to  the  first  and  third  natural  acoustic 
frequencies  of  the  set  up.  The  phenomenon  of 
flamelet  shedding  has  been  observed  at  these 
frequencies.  The  radiation  emission  from  the 
flame,  and  presumably  the  fluctuating  heat  release 
rate,  undergo  phase  changes  along  the  flame  on  the 
scale  of  a  wavelength,  A,  based  on  the  convection 
velocity  of  the  flamelets;  i.e.  A  ~  u  /f  where  uc 
is  the  convection  velocity.  The  acoustic  pres¬ 
sure,  on  the  other  hand,  undergoes  phase  changes 
based  on  the  speed  of  sound.  Thus,  in  certain 
regions  of  the  flame,  the  heat  release  and  the 
pressure  are  in  phase  and  in  other  regions  they 
are  out  of  phase.  For  the  flamelet  shedding  at  a 
preferred  frequency,  the  amplitude  of  the  radia¬ 
tion  maximizes  in  a  region  where  there  is  a 
component  of  it  in  phase  with  the  pressure,  thus 
satisfying  Rayleigh's  criterion  for  flame  driving 
of  the  acoustics.  On  the  other  hand,  the  flame  in 
the  tube  with  the  shorter  exhaust,  is  able  to 
respond  to  the  imposed  excitation  at  the  driving 
frequency.  The  radiation  signals  from  the  flame 
are  higher  when  the  stabilizing  wire  is  near  a 
pressure  minimum. 

Consideration  «f  the  acoustic  velocity  and 
its  relation  to  the  oressure  and  heat  release  is 
also  believed  to  be  important.  Some  of  the 
limitations  of  current  theoretical  models  were 
also  underlined. 
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Fig.  2  Shadowgraph  Film  Sequence 
Showing  Flamelet  Shedding  at  a 
Preferred  Frequency.  Note  Non- 
Negligible  Thickness  and  Structure 
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Fig.  3  Comparison  of  Emitted  CH  Radiation  Amplitudes 
of  First  Natural  Frequency  and  at  Driving  Frequency 
of  385  Hz.  The  Four  Cases  correspond  to  the  Wire  at 
a  Maximum,  a  Minimum,  Between  a  Minimum  and 
Maximum  and  Between  a  Maximum  and  Minimum  for  the 
Standing  Wove  at  385  Hz. 
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Fig.  6  Axial  Variation  of  Phase  of  CH 
Radiation  at  first  Natural  Frequency 
under  Acoustic  Excitation  at  385  Hz. 
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ABSTRACT 

Combustion  instability  problems  often  occur  in 
coaxial  dump  type  ramjet  combustors.  A  potentially 
severe  type  of  instability  is  characterized  by  low 
frequency  (in  the  range  up  to  500  Hz)  pressure  and 
velocity  oscillations.  This  paper  is  concerned  with 
this  type  of  instability.  It  describes  an  investi¬ 
gation  of  the  coupling  between  cone  shaped  flames, 
often  encountered  in  dump  type  ramjet  combustors 
and  longitudinal  acoustic  fields  associated  with 
the  low  frequency  instability.  Specifically  the 
unsteady  behavior  of  a  V-shaped  flame,  stabilized 
in  a  rectangular  combustor,  in  various  longitudinal 
acoustic  fields  is  studied.  A  variety  of  experi¬ 
mental  techniques,  including  high  speed  Schtieren 
and  shadow  photography,  CH  radiation  emission  and 
acoustic  pressure  measurements  have  been  employed 
in  this  study.  It  is  found  that  the  shear  layer 
formed  by  the  flame  exhibits  a  marked  fluid 
dynamical  instability  at  low  frequencies  resulting 
in  the  formation  of  burning  vortical  structures 
that  are  convected  along  the  flame  front.  The 
resulting  unsteady  combustion  process  acts  as  a 
strong  acoustic  source  resulting  in  longitudinal 
pressure  oscillations  This  problem,  therefore, 
presents  an  intricate  coupling  among  the  three 
basic  modes  of  fluid  motion;  that  is.  the  vortical, 
acoustic  and  thermal  modes  A  simple  model  develop¬ 
ed  to  predict  flame  generated  acoustic  spectra  from 
measurements  of  the  unsteady  heat  release  rate 
yields  the  coupling  between  the  thermal  and 
acoustic  modes  The  unsteady  vortical  behavior  of 
the  flame  is.  however,  extremely  complicated  and 
exhibits  many  nonlinear  features,  'n  particular  the 
presence  of  subharmonics  and  harmonics  of  the  main 
response  frequency 

INTRODUCTION 

Recent  developments  of  coax'a'  dump  type 
ramjet  engines  have  been  hindered  by  occurrences  of 
destructive  combust’on  instabilities  '  A  poten¬ 
tially  severe  type  of  instability  's  character i zed 
by  low  frequency  (in  the  range  of  up  to  500  Hz  | 
pressure  and  velocity  oscillations  rhis  type  of 
instability  can  interfere  with  the  shock  system  in 
tho  ramjet  inlet  resulting  in  loss  of  engine  per- 
formance  as  well  as  imposing  excessive  vibrational 
loads  on  the  system.  The  present  paper  investigates 
this  low  frequency  type  of  combustion  instability 

In  general,  the  unsteady  combustion  process 
provides  the  energy  required  for  the  excitation  and 
maintenance  of  the  observed  oscillations  According 
to  Rayleigh's  criterion  ,  driving  of  the  pressure 
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oscillations  by  the  combustion  process  occurs  if 
the  oscillating  heat  release  from  the  combustion 
process  is  in  phase  with  the  pressure  oscillations. 
From  this  viewpoint,  it  becomes  essentia!  to  deter¬ 
mine  the  reasons  and  the  events  that  conspire  in  a 
particular  combustion  instability  problem  to  create 
an  environment  where  Rayleigh's  criterion  is 
satisfied. 

To  date,  a  limited  number  of  experimental 
studies  of  the  coupling  between  the  combustion 
processes  and  acoustic  oscillations  in  a  ramjet 
engine  like  environment  have  been jconducted.  Recent 
studies  include  those  of  Davis  and  Heitor  et 
ar  .  Davis  used  high  speed  photography  to 
identify  two  modes  of  low  frequency  instability  in 
dump  type  combustors.  In  the  first  mode,  the  entire 
combustion  zone  underwent  a  cyclic  oscillation 
while  in  the  second  mode  regular  shedding  of  hot 
spots  from  the  recirculation  zone  at  the  dump  plane 
was  observed.  Heitor  et  al  measured  the  frequency 
and  strength  of  combustion  induced  oscillations  for 
premixed  methane  air  flames  stabilized  on  baffles 
located  on  the  axis  of  a  closed-open  pipe.  They 
found  that  instability  involving  the  fundamental 
acoustic  mode  (i  e.,  a  quarter  wave  type)  occurred 
for  a  wide  range  of  fuel  ip  air  ratios  More 
recently,  Smith  and  Zukoski  studied  combustion 
instability  involving  a  flame  stabilized  behind  a 
rearward  facing  step  They  report  that  under 
conditions  of  instability  large  vortical  structures 
were  formed  downstream  of  the  step  at  an  acoustic 
resonant  frequency  These  were  subsequently  convec¬ 
ted  away  from  the  step  and  the  resulting  unsteady 
combustion  was  found  to  feed  energy  ’nto  the 
acoustic  oscillations  Kel'er  et  a'  idenf‘'ed 
three  mooes  of  instability  in  a  si--lar  set  up  "i 
an  effort  to  understand  the  mechanisms  responsible 
for  flashback  They  attributed  this  phenomenon  to 
the  action  of  vortices  in  the  rec I rcul at 'on  zone  of 
the  step.  These  studies  indicate  that  the  fiu’d 
mechanics  of  the  flame  stabi 1 izat ion  region  plays 
an  important  role  in  combustion  instability 

The  present  study  also  confirms  the  importance 
of  the  fluid  mechanics  of  the  fiax*e  region  The  ex¬ 
periments  earned  out  and  described  later  indicate 
the  presence  of  a  marked  fluid  dynamical  instabil¬ 
ity  of  the  f'ame  region  at  'ow  frequencies  This 
results  m  the  formation  of  burning  vortices  that 
are  shed  from  the  flame  holding  region  and  are 
convected  a'ong  the  flame  front  The  resulting 
unsteady  combustion  acts  as  a  strong  acoust’c 
source  resulting  in  longitudinal  pressure  osc'i'a- 
tions  The  combustion  instability  prob'em  'n  r*m~ 
jets  therefore  involves  an  intricate  coupling  among 
the  three  basic  xtodes  of  f'u'd  motion  that  's  th* 
vortical  acoustic  and  thermal  modes 

EIPERIMENTAL  set  up 

The  set  up  described  m  fig  I  has  been 
developed  It  consists  of  an  inlet,  combustor  and 
exhaust  sections  The  inlet  contains  a  movable 
injector  through  which  a  mixture  of  propane  and  a'r 


is  introduced  into  the  combustor.  The  flame  is 
stabilized  in  the  7.5  x  5  cnr  combustor  section  on 
a  0.8  mm  diameter  nichrome  wire.  The  wire  is 
attached  to  the  combustor  windows  at  half  combustor 
height  and  it  is  heated  electrically  to  improve  the 
steadiness  of  the  stabilized  V-shaped  flame.  The 
exhaust  section  is  equipped  with  two  acoustic 
drivers  which  are  used  to  excite  a  standing 
acoustic  wave  of  desired  amplitude  and  frequency  in 
the  set  up.  The  movable  injector  provides  a  capa¬ 
bility  for  “placing"  the  stabilizing  wire  on  any 
part  of  the  excited  standing  wave;  that  is,  at  a 
pressure  maximum,  minimum  or  in  between.  The  maxi¬ 
mum  length,  L,  available  between  the  injector  face 
and  the  exhaust  section  exit  plane  is  3  meters.  The 
flow  approaching  the  stabilizing  wire  is  parallel 
and  uniform  and  disturbances  are  damped  out  by 
means  of  a  fine  wire  mesh  grid  located  8.  cm. 
upstream  of  the  stabilizing  wire.  The  cold  flow 
Reynolds  number  is  kept  below  10,000  in  the 
Interest  of  maintaining  a  disturbance  free  flow. 
The  combustor  walls  are  water  cooled  enabling  wall 
mounted  pressure  transducers  to  be  used  to  monitor 
the  acoustic  pressure  field.  Capabilities  for  time 
and  space  resolved  measurements  of  CH  species 
radiation  from  the  flame  have  also  been  developed. 
The  concentrations  at, these  species  are  a  measure 
of  the  reaction  rate  and,  thus,  the  heat  release 
rate  and  are  useful  in  describing  the  unsteady 
combustion  field. 

The  acquired  pressure  and  radiation  signals 
are  digitized  and  stored  in  computer  memory  prior 
to  their  Fourier  analyses.  The  developed  Fourier 
analysis  program  yields  both  auto  spectra  and  cross 
spectra  of  the  signals.  This  facilitates  determina¬ 
tion  of  relationships  between  the  pressure  and  the 
unsteady  heat  release  including  their  r«lative 
pnases.  Also,  capabilities  for  identifying  specific 
frequency  components  and  their  amplitude  and  phase 
behavior  have  also  been  developed. 


EXPERIMENTAL  RESULTS 

The  initial  experiments  performed  including 
flow  visual ization  (described  in  Ref.  9)  nad  sug¬ 
gested  a  vortical  instability  of  the  flame  charac¬ 
terized  by  the  Shedding  of  burning  vortices  fro* 
the  flame  holding  wire  around  the  first  natural 
frequency,  f  .  of  the  set  uo  (ranging  between  55-90 
Hz  for  thf  different  effective  tube  lengths 
employed)  A  representative  frame  fro*  the  high 
speed  (6000  frames/ second)  shadowgraph  movie  ’s 
reproduced  m  fig  <?  The  non-negl  tg'bie  thickness 
and  vortex  Mke  structure  of  the  f'eme  region 
should  be  noted  To  classify  this  aspect.  Cm 
rad’et'on  emitted  by  the  fame  (re oresentatwe  of 
the  unsteady  neat  re'easei  was  measured  at 
different  a«-a!  'ocat'ons  along  the  fame  front 

Before  describing  the  experiments.  ’t  s 
appropriate  to  consider  some  features  o'  shear 
’ayer  'nstab'i'ty  Th>$  is  a  complex  orocess,  even 
m  ' ncompress ’ bl e  '’ows  where  a  vasl,  .amount  o'  data 
mey  be  found  m  the  'terature  these  data 

indicate  the  ore%e«ce  of  a  'avored  frequency  range 
'Strouha  number  range)  where  instability  occurs 
Whereas  the  primary  energy  causing  f’ow  instability 
com**  'rom  the  mean  flow  v'a  interaction  w‘th  the 
Reynolds  stresses,  nominee'  ’nterect’ons.  e  g 
merging  and  breakuo  of  the  'ormed  vo't'ces  a'so 
cause  energy  transfer  among  different  frequency 
components  This  results  •«  the  generation  o' 


harmonics  and  subharmonics  of  the  primary 
instability  frequency. 

For  the  flame  experiments,  the  preliminary 
investigation  had  shown  that  the  preferred 
frequency  for  the  flamelet  s., adding  was  around  f  . 
Hence,  the  Fourier  analysis  program  for  the 
radiation  data  was  extended  to  pick  out  the 
frequency  components  not  only  at  f  but  also  at 


a) 

V2 

the 

first  subharmonic 

O 

o 

b) 

2  fo 
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first  harmonic  of 

f  and 

0 

c) 

fd 

the 

external  driving 

frequency,  if 

any. 


A  typical  result  is  plotted  in  Fig.  3  which 
shows  the  dependence  of  the  radiation  amplitude 
upon  the  axial  distance  x  from  the  flame  holding 
wire.  In  this  case  there  was  no  external  driving 
and  f  “  80  Hz.  Note  the  sharp  increase  in 
amplitOde  at  all  three  frequencies  ( i . e .  ,  80,  40  i 
160  Hz)  just  downstream  of  the  wire.  Whereas  the 
signal  at  80  Hz  is  dominant,  there  is  a  strong 

signal  also  at  40  Hz.  The  increase  in  amplitude 

with  progressing  axial  distance  that  is  observed  is 
characteristic  of  shear  layer  instability.  The 
signals  start  decaying  towards  the  end  of  the  flame 
(i.e.,  x*20  cm)  which  is  probably  due  to  the  wall 
quenching  effect  on  the  flame  since  the  flame  meets 
the  duct  walls  at  its  downstream  end. 

Figure  4  plots  the  phase  of  the  radiation 
signal  at  80  Hz  with  respect  to  the  acoustic 

pressure  as  a  function  of  x.  Note  the  region  of 
phase  reversal  in  a  small  region  just  downstream  of 
the  wire.  Since  the  phase  curve  slope  is  related  to 
the  convection  velocity  of  the  disturbance,  the 
region  of  phase  reversal  indicates  a  negative 

convection  velocity  and  hence  represents  the 
recirculation  region  downstream  of  the  wire 
Downstream  of  this  region  the  phase  increases 
smoothly  and  monoton i ca 1 1 y .  The  behavior  of  the 

phase  distributions  of  the  other  'ow  frequency 
component  signals  is  similar.  It  is  found  jlso  that 
the  degree  of  phase  variation  over  the  fame  ’s 
larger  for  the  higher  frequency  components  ('  and 
above)  as  may  be  expected,  due  to  their  smj''er 
vortical  wavelengths 

The  convect'on  velocity  m  the  flame  rn'on 
for  the  rjd'ation  signal  at  f  •$  plotted  m  r-g 
5  Note  the  'ncrease  m  the  velocity  w'th  the  a«-a’ 
distance  due  to  the  acce'*rat’on  caused  by  t*e 

density  gradients  The  rat'o  of  about  4  0  ‘o'  the 
convection  velocity  at  the  downstream  end  to  that 
at  the  upstream  end  's  equa'  to  the  temperature 
rat'o  across  the  flame  A'so,  the  va'ue  of  the 

convection  veioc’ty  (1  3  m/sec)  at  the  upstream  end 
is  c'ose  to  the  co'd  f'ow  ve'oc'ty  o'  ;  5  m,  set 
This  suggests  that  the  vort’cai  structures  are 

oas’ca'y  convected  at  the  average  speeds  o'  the 

not  and  co'd  f ’ ows  at  any  st'eamw’se  ’ocat'o",  . 
i'ong  the  fame 

The  rfd'at'on  spectrum  'rom  the  ent're  *  am* 
was  a'so  measured  A  typ'ca'  resu  t  s  shown  *n 
c'g  6  it  *s  seen  to  be  dominated  by  'he  ow 
'requency  components  under  '  Tnus.  at  oca 
stat’ons  a : ong  the  f ame  ''8nt,  'he  rad’afo" 
s'gn*ls  at  frequencies  c'ose  to  '  dominate  but 

considering  the  radiet'on  'rom  the  e8f 'e  '’em*  the 

ow  frequency  components  dominate  Tne  -eason  'o- 


this  is  that  the  radiation  signals  from  different 
parts  of  tha  fla**a  cancal  out  for  tha  higher 
fraquancy  components  dua  to  tha  largar  phasa 
variation  of  tha  signal  ovar  tha  flame.  for  tha  low 
fraquancy  components,  tha  phasa  variation  Is  small 
so  that  cancallatlon  affacts  ara  also  small. 

Tha  affact  of  external  acoustic  excitation  at 
a  frequency  f  may  oe  categorized  into  the  follow¬ 
ing  classes: 

a)  f  =  nf  Where  n  is  an  integer  >  1.  This 

corresponds  to  driving  at  a  harmonic  of  f  . 

B )  f  =  f  /n.  This  corresponds  to  driving  at  a 

subharn&nic  of  f  and 
o 

C)  f  is  arbitrary.  This  Corresponds  to  driving 

at  neither  a  harmonic  nor  a  subharmomc  of  f 

o 

To  investigate  Type  (a)  driving,  f  was  sat  at 
160  Hz  To  investigate  Type  (c)°  different 
frequences  ware  choi*n;  e  g.,  385  Hz  and  270  Hz 
In  these  cases,  tests  ware  run  so  that  the  f'ame 
holamg  wire  was  on  afferent  locations  of  the  ex¬ 
cited  standing  wave.  '  e  .  at  a  pressure  maximum,  a 
pressure  minimum  and  m  Between  unfortunate  1 y . 

Type  (bi  cou'd  not  Be  run  on  the  sat  up  dua  to  the 
inefficiency  o'  the  acoustic  drivers  at  ow 
f requenc i as 


effect  of  one  frequency  component  of  the  radiation 
signal  upon  the  other.  Simitar  types  of  Behavior 
have  bean  reported  in  studies,  of  coherent 
structures  in  nonreacting  flows  .  This  nonlin¬ 
earity  has  been  attributed  to  vortex  merging  and 
pairing  in  the  shear  layer. 

!t  is  clear  that  in  the  reported  experiments 
the  flame  is  vortical ly  unstable  at  the  first 
natural  frequency,  f  ,  of  the  set  up.  To  maintain 
this  instability  a  feedback,  mechanism  through  which 
disturbances  at  f  are  “felt"  By  the  'lame  must  be 
present  in  the  experimental  set  up.  Since  f  is  the 
first  natural  acoustic  frequency  of  the  sA  up  it 
is  natural  to  seek  an  appropriate  acoustic  source 
in  the  set  up  instrumental  in  closing  the  feedback 
loop  between  the  vortical  instability  and  the  duct 
acoustics. 

According  to  Rayleigh’s  criterion.  the 
unsteady  heat  release  from  the  'lame  should  ce  m 
phase  with  the  acoustic  oressure  on  the  average 
over  the  entire  flame  if  it  is  to  drive  the 
oressure  oscillations.  This  criterion  may  Be  stated 
m  the  following  form  in  the  frequency  doma’r 

/  ,,  v,,,. 

i  ,  acic 


for  tests  w’tn  f  of  e’tner  'ype  (ai  or  ry oe 
t  c  '  a  very  interesting  ohenomenon  was  observed 
figure  '  o'ots  the  'ad'at-on  amplitudes  vs  «  at  80 
hz  and  160  hz  when  f  =  160  hz  Note  the  sharp  - ■ se 
’h  the  signal  at  1%0  hz  ust  downstream  of  the 
wire  In  fact  for  i  me  s'gnai  at  160  hz 


dominates  t"e  sgna'  at  80  hz 


However . 

oast  this  axia1  location,  the  80  hzu  i  -  e  .  * 
s ' gna i  aga’n  dominates  the  160  ny  s’gna'  »  s’«’'ir 
behavor  was  observed  -  tests  wt"  ’ype  ic 


frequences  it  was  *ound  that  '"dependent  of  the 


ocat  o"  of  the  *  ame  n0  o-"g  w -e  *"e  e<cted 
acoust’t  wave,  t't  -ad’afo"  s  gna  ,ery  •  ose  to 
the  wre  was  dominated  o*  the  d-'v'ng  '-eque-ry  ‘ 
But  oast  t "  s  '"'ta  'eg'or  *ne  s-gna  at 

dominated  for  «>amp  e  -g  H  lest-'Oes  t"# 

rad'at'on  amo'  tudes  at  80  ny  inj  igs  wy  -us' 
downstream  of  the  x’re  the  s'gna  at  '  dominates 
the  size  of  t"e  'eg'or  dowrst'eam  of  t"j  w  -e  where 
the  S’gna'  at  *  dominates  the  s’gna  at  *  depends 
on  the  amp  ’tude  of  the  d'’v'ng  nrreei^ng  »i>.' 
’nc-ees’"g  amc'Hude  of  dr'v'ng  jng  nr  *_  tse  ' 
idecreasmg  wt"  increasing  * 


0 


thus.  the  ‘ame  may  Be  thought  o*  as 
comprising  two  -eg'ons  e  a  a  -ea-  wa«e 
reg'o«  'espons've  to  Htamj  eif'fet’on  an«  B  a 
f»r  wake  reg'on  where  the  do«'"ant  *-a 


'*Oue"C  y 


-he  effect  o'  utern  ir.'ng 
of  the  radtjt’on  s’gna  at 


Be"av'o' 


'"vest ' gated 


'  gur  es 


o  and  ;c  ^lesc''be  t"e 

amplitude  end  ohase  o'  ’**e  ,  gna  at  80  "t  ' 

with  and  w'thout  dr'v’ng  at  i6C  ny  'Xfry  does  "£t 
appear  to  be  my  s'gn'f'cant  effects  o'  t"e  dr’v'ng 
on  the  phase  behev'or  at  9C  «y  implying  »h$t  •"» 
c  onvec  t  *  on  ve'ocity  o'  th#  d'sturhence  's  naf 
e'tered  But  the  samp  ■  a""ot  be  sa'd  about  t"e 
amp  *  1 1  ude  vif  *  a  t  ’  on  ph  i  1  e  'he  t'a'ti  o'  t  "e  s'gna  s 
>n  both  rases  s  S’m>'e'  '»e  eve  s  •'» 

s*gn<f'cant'y  d’f'erent  oe'ng  sgrh  owe'  ■"  ?"* 

dr  * v i ng  cast  this  emphasizes  the  "on1 inter  nature 
of  the  f'amp'et  t''Odd*ng  phenomenon  'hat  ■%  the 


where  S  *s  the  cross  spectrum  Between  t't 
Bressure°lnd  "eat  re  ease  '  e  rad’at'on)  and 

DO 

■s  the  B"ase  Between  the  unsteady  neat  re'ease  a"d 
the  pnessune  1'  the  cnter-on  s  sacs'-ea.  '"e 
feedback  ’  ocb  Between  the  shear  aye'  -nstaB'  -ty 

o*  the  f'ame  and  t"e  duct  acoust'cs  may  oe  c’osea 

"he  mtegrand  o'  c.o  i  ’s  oott*<d  ■"  f'g  :: 

as  a  'unct'on  o'  i  'or  '  :  SO  i  "o  dr', mg 

eiBer-me"!  ■"»  ,a  u»  o'0  t"e  "ter>  ’s  ‘ouna  tc 
Be  oos'f.o  "  accordance  wth  '."e  '';e''on  4o' 

ot"*r  ow  ''eouency  .o"iponents  se  ow  •  '"e 

r • t e ' ■ o "  s  a  so  ’ c  Be  sat  s‘  ed  '"d'fat  ng 

'"a*  me  '  ame  a  sc  d’-.es  'he  P'tssj'f 
:■  s  ■  at  ’  o"  s  a’  '."ese  frequencies 

■he  sat's'ymg  o'  ’.he  criterion  "owfii'  does 
*ot  e*c  ude  ’"e  o'esenre  o'  ot"er  'mpc>'’ant 

ecwst'c  soj"»s  the  set  up  for  name  » 

mp’ngpmpnt  o'  'he  vont'ces  aga'nst  a  so'-d  su'*are 
e  g  'he  due'  m  's  cou  d  a'so  ontr'bute  to  the 
gene-et’on  'f  t"e  duct  acoustics  as  "as  Bee" 
suggested  to  y  cry  •"  segmented  so  d  orope  an* 

-oc«et  mote's  '"us  •'  Becomes  "*c»ssar  i  to 

B'td'rt  the  acoust'c  C'essu'e  sont'jn  By 
o"s'dermg  'he  unsteady  "eat  **  ease  as  f"e 

pn-mpry  acOuSt’-  source  and  ompj'  ng  f  h 1 s  P'ed'I 
'■on  w't"  the  eipe'impnta  o'essu'e  spect'a  I' 
set  *  s'ac  tory  ag'tement  s  obta‘"ed  t"en  t"e 

<"steady  "ea’  -e  ease  -ate  may  Be  'n»"  to  he  t"e 
on  y  moo'tant  acoust'i  sourre  'h-s  o'ed'ct-on  -s 
desc''Bed  ■"  t»e  »eit  section  and  't  shows  that 

'h-s  "deed  appea's  to  Be  'he  aye 

’“f  £>•{  ’  ICAi  CCMKSIDf  »At  !0*(S 
»« per ’went |  -esu  ts  lest ’’bed  m  to# 

O'freed'ng  sect'on  demonstrate  that  ombust’o" 
"stab"  ’y  ■ "  me  set  cp  s  onnec  ted  wt"  'he 

v  o '  t  ’  c »  msteb  'ty  o'  the  *  ’  amp  'he  density  jump 

ac'ots  thp  'eme  auses  sha'p  veoc'ty  g'ad'ents 


resulting  In  an  intansa  shear  layer.  This  shear 
layer  is  vortically  unstable  in  a  certain  frequency 
range  depending  upon  the  flow  conditions.  For  the 
experimental  conditions  considered  here,  this  range 
appears  to  be  limited  to  frequencies  below  100  Hz. 
Unsteady  combustion  in  the  vortices  formed  due  to 
the  shear  layer  instability  can,  in  accordance  with 
Rayleigh's  criterion,  generate  an  acoustic  field. 
If  the  shear  layer  instability  frequency  range 
contains  a  natural  acoustic  frequency  of  the 
system,  then  the  entire  instability  process  is 
reinforced  and  may  result  in  large  pressure 
fluctuation  levels. 

A  comprehensive  theoretical  analysis  of  this 
problem  is  still  beyond  the  state  of  the  art 
modelling  techniques.  Since  the  vortical  structures 
appear  to  be  Ouasi  deterministic,  traditional 
statistical  analyses  which  assume  randomness  cannot 
reveal  completely  all  the  features  of  the  flow. 
Recent  trends  have  been  towards  a  direct  numerical 
simulation  .Jif  vortical  structures  in  mixing 
layers.  '  '  However,  these  methods  are 
restricted  to  low  or  moderate  Reynolds  number 
because  of  the  limited  range  of  temporal  and 
spatial  scales  that  can  be  resolved  with  today's 
computers.  As  a  result,  only  the  global  features  of 
the  flame  may  be  reasonably  modelled. 

As  a  first  step  in  this  direction  a  simple 
model  is  developed  to  predict  the  pressure  spectrum 
from  the  unsteady  heat  release  rate  spectrum.  If 
satisfactory  agreement  is  obtained  it  would  also 
indicate  that  the  unsteady  heat  release  rate  is  the 
primary  acoustic  source  in  the  set  up  as  mentioned 
in  the  preceding  section.  In  the  analysis  the 
unsteady  heat  release  rate  is  specified  by  the 
unsteady  CH  radiation  emitted  by  the  flame 
(obtained  experimentally).  A  future  goal  would  be 
to  oredict  the  unsteady  heat  release  rate  also, 
oemaps  from  knowledge  of  the  steady  flow 
conditions 

The  analysis  considers  only  low  frequencies 
(up  to  300  Hz  approximately].  In  tils  range.  the 
flame  zone  may  be  assumed  compact  compared  to  the 
acoustic  wavelength  Thus,  the  flame  is  assumed  to 
be  localized  at  x  =  I .  (  i th  x  =  0  being  the 
injector  and  i  i  l  the  exhaust  plane)  The  steady 
state  temperature  ’s  assumed  to  be  T  in  the  region 
0-  »  •  i ,  and  t,  (burned  gas  temo^v-ature )  m  the 
region  l  4>  ■  •  \  T  is  taken  to  be  an  average 
downstrew  temperature^  as  there  are  some  neat 
'osses  from  the  water  coo'ed  walls  of  the  duct 

As  only  plane  acoustic  waves  are  of  interest, 
only  the  one-d  imens '  ona  i  wave  equat'on  for  the 
acoust'c  pressure  p  s  considered  moreover, 
steady  flow  ve'oc’ty  effects  are  neglected  due  to 
the  extremely  'ow  Wech  number  of  the  f'ow  i0  005 
approx  imate  ’  y  I  for  anal  temperature  grao'ents 
w<th  unsteady  neat  release  q'(x)  this  wave  equation 
may  be  wntten  m  the  frequency  domain  as 


Gix.c  l  ’s  obtained  from  Ea  (3)  by  sett'ng  y  r  L 
It  ma£  oe  constructed  by  standard  methods  '  ’he 
comparison  with  the  expenmenta'  pressure  spectrum 
x i i '  be  made  at  the  injector  'he  re'evant  Green's 
function  Gio.tji  turns  out  to  be 


where  i  i  x,  „ .  and  s  - 

A  A 

anq  ‘  k i  *  and  c ,  ana  c  , 

t  * 

are  the  average  speeds  of  sound  upstream  ana 
downstream  o'  the  f’ame  respe r t ’ ve 5 y 


Th«  pressure  spectrum  S  is  obtained  by 
multiplying  Eq.  (4)  by  its  coR&lex  conjugate  and 
talcing  an  ensemble  average.  The  result  is 


a  ,  igco.i.^:2  s^u  =  Uj_) 


where  S  is  the  net  unsteady  heat  release  spectrum 
from  th8qentire  flame.  It  is  further  assumed  that 


where  $  is  the  CH  radiation  autospectrum  measured 
from  the  entire  flame  and  B  is  a  positive  real 

constant  relating  the  magnitudes  of  the  two 
spectra. 

A  typical  comparison  of  the  predicted  and 

measured  pressure  spectra  is  shown  in  Fig.  12.  The 
agreement  is  good.  It  is  also  seen  that  while  the 
spectrum  shows  a  peak  at  f  (80  Hz)  as  expected, 
there  are  also  substantial  pressure  levels  at  lower 
frequencies  indicating  that  the  flame  also  drives 
at  these  lower  frequencies.  As  noted  in  the 

previous  section,  Rayleigh's  criterion  is  also 

satisfied  at  these  lower  frequencies. 

It  should  be  noted  that  the  magnitude  of  the 
denominator  of  the  Green's  function  G(Eq.  5)  is 
minimized  for  the  natural  acoustic  frequencies  of 
the  duct.  This  means  that  driving  of  the  natural 
acoustic  frequencies  is  favored  as  is  expected 
Thus,  if  high  levels  of  unsteady  heat  release  rate 
are  obtained  at  a  natural  frequency  of  the  duct, 
extremely  large  pressure  oscillations  may  result  as 
is  evident  from  Eq.  (6) . 

In  the  experiments  described  herein,  it  is 
be'ieved  that  the  fundamental  acoustic  frequency  f 
of  the  set  up  also  belonged  in  the  range  0? 
unstable  frequencies  for  the  flame  shear  layer 
instability  An  interesting  question  that  may  be 
raised  is  whether  the  acoustic  natural  frequencies 
influence  m  some  way  the  unstable  frequency  range 
for  the  shear  layer  instability  by  shifting,  for 
example,  this  frequency  range  To  answer  these 
types  of  questions  it  is  o'anned  to  repeat  some  of 
the  described  tests  m  a  shorter  version  of  the 
developed  set  up  This  will  serve  to  increase  the 
first  natural  acoustic  frequency  and.  presumably, 
remove  it  from  the  unstable  frequency  range  for  the 
flame  shear  layer  instability 

CONCLUSION 

The  response  of  wedge  shaoed  flames,  similar 
to  those  often  occurring  m  ramjet  combustors,  to 
acoustic  excitation  nas  been  studied  The 
investigation  has  shown  that  the  unsteady 
combustion  in  vortices  shed  from  the  flame  holding 
region  is  the  major  cause  of  the  observed  pressure 
oscillations  m  the  set  up,  and  could  also  serve  as 
the  driving  mechanism  for  low  frequency  ramjet 
instabilities  The  vortical  structures  arise  due  to 
the  instability  of  the  shear  layer  formed  by  the 
flame  This  instability  and  therefore  the  vortical 
structures  are  restricted  to  a  limited  frequency 
range  Thus,  significant  rates  of  unsteady  heat 
release  (and  driving)  are  also  confined  to  the  same 
frequency  range  Overlap  of  an  acoustic  resonant 
frequency  with  the  frequency  range  of  the  flame 


shear  layer  instability  can  result  in  extremely 
large  pressure  oscillations  and  should  be  avoided. 
Better  understanding  of  the  origin  and  detailed 
structure  of  the  vortices  is  the  key  to  their 
control  and,  perhaps,  the  elimination  of  combustion 
instability  in  ramjets. 
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.  SOUND  GENERATION  BY  DUCTED  FLAMES 


*  **  *** 
U.G.  Hegde  ,  D.  Reuter  and  B.T.  Zinn 

Georgia  Institute  of  Technology 

Atlanta,  Georgia  30332 


ABSTRACT 


The  sound  field  established  by  a  v-shaped  flame  confined  in  a 
rectangular  duct  is  investigated  both  theoretically  and  experimentally.  A 
theoretical  model  is  developed  to  predict  pressure  spectra  caused  by  the 
unsteady  heat  release  from  the  flame.  Comparisons  between  the  theoretical  and 
experimentally  measured  spectra  confirm  the  validity  of  £he  model  and  are 
also  reported.  It  is  also  found  that  changes  in  the  flow  field  in  the  flame 
zone  can  significantly  modify  combustion  rates.  These  modifications,  in  turn, 
affect  the  generated  sound  field  and  are  important  in  determining  the 
pressure  levels  in  the  duct.  The  results  of  this  investigation  are  applicable 
to  combustion  noise  and  instability  studies  in  a  variety  of  burner  and 
propulsion  system  configurations. 
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NOMENCLATURE 

duct  width 
duct  height 

constant  defined  by  Eq.  (17) 
constant  defined  by  Eq.  (18) 
speed  of  sound 

specific  heat  at  constant  pressure 
Green's  function 
modified  wave  number 
length  of  duct 

length  of  duct  from  entrance  to  flame  2one 

steady  state  Mach  number  at  duct  inlet 

steady  state  Mach  number  at  duct  exit 

heat  release  rate  per  unit  volume 

gas  constant 

pressure  autospectrum 

unsteady  heat  release  rate  autospectrum 

CH  radiation  autospectrum 

time 

temperature 


steady  state  cold  gas  temperature 
steady  state  f ’ amt  temperature 
velocity 

a  < i a  1  coord'nate 

aula1  sepa  at’on  of  f ' ame  holders 
jtrT-ra'  roord,riate 

/e'^'a’  separaT'on  o*  * 1  am*  *‘</d*«,s 


y.'s.  v-  v.  v  ■» 


2 


transverse  coordinate 


2  ,2  unmodified  and  modified  specific  acoustic  impedances  at  entrance  to 
o  o 

duct 

A 

2^,2l  unmodified  and  modified  specific  acoustic  impedances  at  exit  of 
duct 

a  defined  by  Eq.  (14) 

0  defined  by  Eq.  ( 14) 

specific  acoustic  admittance  of  duct  side  walls 
6  Dirac  delta  function 

y  ratio  of  specific  heats 

uj  angular  frequency 

X  wall  loss  factor  (per  unit  length) 

x  acoustic  wavelength 

Superscripts 

time  or  ensemble  average 
1  time  dependent  part  (fluctuation) 

*  complex  conjugate 


Subscripts 

^  Fourier  component 

1  relating  to  the  region  upstream  of  the  combustion  zone 

2  relating  to  the  region  downstream  of  the  combustion  zon* 
vector  quantity 

Other 


cross  sect'on  average 


introduction 


Sound  generated  by  unsteady  combustion  orocesses  is  of  interest  in 

several  applications  including  propulsion  devices,  industrial  furnaces, 

burners,  heaters  and  so  on.  In  the  majority  of  cases,  the  generated  sound  is 

undesirable,  as  for  example,  in  ramjet  engines  where  the  problem  may  become 

so  severe  that  unstable  operation  of  the  engine  may  occur*.  In  some  cases, 

however,  as  in  pulse  combustors,  sound  pressure  variations  are  essential  to 

2 

maintain  operation  . 

Unsteady  combustion  processes  generate  sound  as  a  result  of  unsteady 
expansion  of  the  fluid  undergoing  reaction.  Depending  upon  the  nature  of  the 
f 1  ow  field  and  the  combustion  process,  unsteady  combustion  may  arise  in 
different  ways.  For  example,  the  injection  of  fuel  or  combustible  mixture  may 
be  unsteady,  the  flow  field  may  be  turbulent  or  there  may  be  feedback  of  the 
acoustic  motions  on  the  flame. 

It  is  convenient  for  the  purposes  of  this  paper  to  categorize  confined 
blames  as  those  burning  in  an  enclosure  and  unconfined  flames  as  those 
burning  m  a  free  space.  The  distinction  is  important  in  characterizing  the 
generated  sound  field  Thus,  for  unconfined  flames,  the  raciated  sound  field 
is.  generally,  non  directional  and  monopole  in  character^.  For  confined 
fames,  on  the  other  hand,  reflections  from  the  boundaries  and  wave  guide 
ef#ects  can  cause  preferred  directionality  of  the  acoustic  motions  as  well  as 
sustain  discrete  frequency  or  narrow  band  excitation  around  the  natural 
frequencies  of  the  enclosure 

There  have  been  several  recent  investigations  of  combustion  generated 
5  4 

sound  Strahle  '  has  developed  a  general  theory  of  combustion  generated 
no'se  for  free  ♦'ames  m  the  low  frequency  M m i t  and  shown  that  the  radiated 


sound  pressure  is  related  to  the  time  derivative  of  the  heat  release  rate 
averaged  over  the  combustion  region.  The  generated  sound  in  such  cases  is 
generally  broadband  with  a  frequency  content  related  to  the  temporal 
variations  in  volume  of  the  combustion  region.  Experimental  verification  of 
the  theory  has  been  obtained  by  Sivasankara  et  al^. 


Sound  generated  by  confined  flames  have  been  studied  in  relation  to 


combustion  instabilities  in  rocket  motors  and  ramjets.  A  review  is  provided 


by  Barrere  and  Williams”.  A  recent  experimental  investigation  of  noise 


sources  in  a  ramjet  like  combustor  is  presented  by  Poinsot  et  ar  .  The 


investigation  presented  here  is  also  part  of  a  larger  research  effort  on 
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combustion  instabilities  in  dump  type  ramjet  combustors  . 

In  this  paper,  a  theory  for  predicting  the  sound  generated  by  flames 
enclosed  in  a  duct  is  developed.  Knowledge  of  the  unsteady  combustion  field 
is  presumed  and  attention  is  focused  on  longitudinal  wave  motions.  This  is 
presented  in  the  Theoretical  Considerations  Section.  The  theory  is  tested  in 
an  experimental  set  up  in  which  a  v-shaped  flame  is  stabilized  in  a 
rectangular  duct.  Measurements  are  made  of  the  unsteady  heat  release  from  the 
flame,  in  terms  of  emitted  radiation,  for  input  into  the  model.  The 
calculated  pressure  spectra  in  the  duct  are  then  compared  with  experimentally 
measured  spectra  to  test  the  theory.  The  results  are  presented  in  the 
Experimental  Efforts  Section.  In  addition,  the  relationship  between  the 
frequency  content  of  the  unsteady  combustion  and  the  generated  sound  is 
investigated.  It  is  also  shown  that  the  combustion  rates  are  sensitive  to 
changes  in  the  flow  field  in  the  flame  region  and,  in  turn,  affect  the 
generated  sound  levels. 


THEORETICAL  CONSIDERATIONS 


Consider  a  rectangular  duct  similar  to  the  one  in  Fig.  1.  The  x-axis  is 
parallel  to  the  duct  axis  and  the  cross  sectional  dimensions  are  a  and  b.  The 
length  of  the  duct  is  L.  A  combustible  mixture  is  introduced  at  x=0  and  the 

combustion  occurs  in  a  region  around  L^.  The  mach  number  of  the  flow  is 

"2 

assumed  to  be  small (i.e.,  M  «1).  A  typical  axial  temperature  profile  which 

O 

might  be  present  in  the  duct  is  also  shown  in  Fig.  1.  Due  to  heat  conduction 
upstream  of  the  combustion  region  there  is  a  gradual  rise  in  temperature  from 
the  cold  temperature  Tq  at  x=0.  A  sharp  rise  in  temperature  occurs  in  the 
combustion  zone  to  the  flame  temperature  T^  .  Downstream  of  the  combustion 
region  the  temperature  may  again  drop  due  to  heat  losses  to  the  duct  walls. 

The  wave  equation  for  the  pressure  fluctuations  in  the  duct  will  be 
briefly  formulated.  The  pressure  fluctuations  will  be  assumed  to  be  small 
with  respect  to  the  mean  (i.e.,  steady  state)  pressure  levels  allowing  the 
application  of  the  linearized  conservation  equations  to  obtain  the  required 
solutions.  In  addition,  the  actual  details  of  the  mean  velocity  field  in  the 
duct  will  be  ignored.  There  are  two  reasons  why  this  is  permissible.  First, 
the  mach  number  is  taken  to  be  small  so  that  convective  effects  in  the  wave 
equation  may  be  ignored.  Secondly,  although  the  flow  field  is  instrumental  in 
determining  the  combustion  rates,  the  actual  details  of  this  process  is  not 
considered  herein  and  the  needed  reaction  rates  will  be  assumed  known.  This 
means  that  it  is  expected  that  the  effects  of  the  flow  field  on  the  generated 
sound  will  be  obtained  implicitly  through  the  unsteady  combustion  process  and 
there  will  be  no  significant  effects  on  the  sound  field  only  by  virtue  of 
changes  in  the  flow  field  which  do  not  affect  the  reaction  rates.  This 
approach  is  justified  by  the  fact  that  unsteady  combustion  is  a  monopole  type 


of  acoustic  source  and  is  a  much  more  efficient  sound  source  than  pure  flow 
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noise  sources  which  are  either  quadropole  or  dipole  in  character  .  In 
addition,  experimental  results  obtained  as  part  of  this  investigation,  and 
discussed  in  the  next  section,  are  in  agreement  with  the  model  predictions. 

Neglecting  viscous  and  heat  conduction  effects  (although  mean 
temperature  and  density  gradients  are  allowed)  and  assuming  a  perfect  gas 
behavior,  the  linearized  conservation  equations  can  be  expressed  in  the 
following  form: 

Continuity: 


If  +  » 


'  +  u'.vj  -  0 


Momentum: 


Energy : 


State: 


°  57  ■  ' 


;cr  £'♦ 


P'  -  RT  o '  +  Re  T' 

Equation  (1)  may  be  rewritten  in  terms  of  p',u'  and  q'  by  substituting 
for  o'  and  T'  from  equations  (4)  and  (3)  and  making  use  of  the  steady  state 


equation  of  state 


p  ■  R  p  T 


yielding 


1  3p/  -  >  g ' 

—y  +  p  V.u  * 


(6) 


WT  w*  m  f*w*w 


where  c2  -  yr* 

Of  Eq  (2)  and  S“l>tract1ns  the  0»rt<4'  t,m* 
Taking  th.  divergence  °f  E<,(Z' 

derivative  of  Eq.(«)  •»*  obt*',1S 


,  ,,  i  iV.-l  &' 

7.(*  CP  )  -  2  .  2  C  at 

C  oC  r 
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Mnn  for  the  pressure  talcing  into 

Equation  ,7)  *  «•  . . .  process  ,„q  the  spatial 

account  the  heat  release  from  the  co  idered  so  that  the 

variations  of  the  temperature.  Harmonic  motion 

fluctuations  may  be  written  in  the  form 

.  .  p  ^  and  q’  -  P, 
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tion  „  now  focussed  only  on  longitudinal  The  axial 

^  ^  •  be  obtained  by  integrating  Eq.(8)  over  the  duct  cross 

wave  equation  may  b  assumed 

•  that  T/c2  is  independent  of  the  location,  it  wi 
section.  Noting  ,on8itudinal  wave  motions,  it 

that  it  may  be  replaced  by  <  T>/<c>.  Also, 

is  permissible  to  replace 
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i.  (ft  ^ 


J..  A-r  hv  T 


d<p  > 
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as 


'  01 


<  p 


U) 


Then,  for  longitudinal  wave  motion  Eg . { 8 )  may  be  expressed  as  *o  ows 


7-<?>j-<p>  +  <  f  >  <k  2  >  <  p  >  -i 
dx  dx  tu  oo  ii 


r  <q  > 

L  HI 


4  ^ 


whene,  <koo2>-  [  -  t«>  ^  ^lir^ 


<c  > 


and  ii>  is  the  specific  acoustic  admittance  of  the  side  walls  of  the  duct  and 
is  assumed  to  be  independent  of  location  along  the  duct.  For  rig;d  wa'led 
ducts,  this  admittance  is  zero.  For  nearly  rigid  walls,  as  is  cons;dered 
here,  the  admittance  remains  small  and  the  modified  wave  number  <kQ0>  may  be 
approximated  by^ 


<  k  > 
oo 


Hi  . 

<  c  > 


i  X 


where  x  Is  known  as  the  wall  loss  factor. 

Removing  the  now  superfluous  <  >,  Eq.  (9)  may  be  written  as 


a  dp 
i-  t 

dx  dx 


+  Tk 


oo 


2  zim 

PHj  =  C 
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(10) 


where  it  is  now  understood  that  cross  sectional  averages  of  quantities  are 
under  consideration. 

The  boundary  conditions  at  x=0  and  x=L  may  be  written  in  terms  of  the 
specific  acoustic  impedances  ZQ  and  ZL  at  these  two  planes.  It  should  be 
noted  that  although  convection  effects  may  be  negligible  in  determining  the 
wave  structure  at  low  mach  numbers,  convection  losses  at  the  end  planes  may 
be  important  in  determining  the  amplitudes  of  the  motions.  This  may  be 
incorporated  as  shown  in  Ref.  11  and  the  relevant  boundary  conditions  are 


along  with  the  same  boundary  cord;tions  as  for  p  .  sc’^t'C"  ‘:r  r  ~av 

U*  •* 

then  be  written  as 


G(x,x  )  dx 
o  o 


(12) 


An  analytical  solution  will  now  be  obtained  for  the  important  case  when 
the  length  of  the  combustion  zone,  AL^ ,  is  much  smaller  than  the  acoustic 
wavelength, A.  In  such  cases,  the  combustion  zone  may  be  treated  as  a 
discontinuity  at  x-Lj  as  far  as  the  acoustics  are  concerned.  The  unsteady 
heat  release  rate  q^(x)  may  then  be  written  as 

qu,  3  Vx)  6  (x  '  V 


and  the  corresponding  solution  for  p  becomes 


(13) 


P  (x) 
jj 


igiq^d^)  C  (X>L| ) 


^e  Green's  funnction  G(x.L^)  may  be  obtained  from  tq.(ll)  by  setting 
<3  =  lj.  ^he  required  jump  conditions  across  the  discontinuity  at  are 


■  G  's  continuous  at  and 


i  1 


x 

dx 


*L- 


The  second  condition  is  obtained  by  integrating  Eq.(ll)  across  L^. 

In  general,  most  of  the  steady  state  temperature  change  will  occur  in 
the  combustion  region  as  shown  in  Fig.  1.  This  change  is  accounted  for  in  the 
matching  conditions  across  the  flame  zone.  The  changes  in  the  steady  state 
temperature  away  from  the  combustion  zone  will  now  be  accommodated  in  an 
approximate  fashion  to  simplify  the  computations.  Upstream  and  downstream  of 
the  combustion  zone,  constant  (although  different)  average  temperatures  will 

s 

be  prescribed.  Thus,  in  the  region  0  <  x  <  l^,  the  steady  state . temperature 
is  taken  to  be  T^  (approximately  equal  to  the  cold  gas  temperature)  and  in 
the  region  <  x  <  L  ,  the  steady  state  temperature  is  taken  to  be  T? 
(approximately  equal  to  T^,  the  flame  temperature,  provided  heat  losses  to 
the  walls  are  not  significant). 

The  comparison  with  the  experimental  pressure  spectrum  will  be  made  at 
x  =0.  The  Green's  function  G(0,L^  )  may  be  obtained  by  standard  methods*'  a-o 


it  turns  out  to  be 


cos  0  +  i  sin 


lZoJlZL 


|k^T^[Zosina  +  i  cos  or]  [ZLcos  0t  i  sin  p]  ‘k2T2^“ZLSin^  +  i  cos  0]  [Zqcos  or  -  i  s 

(14) 

where  a  =  (kQ0)1  4  ,  0  *  (kQ0)2  (4  -  4)  »  4  "  ^oo*!  *  k2  “  ^oo^ 


with  (k  4  and  (kQ0)2  being  the  relevant  wave  numbers  upstream  and 

A  A 

downstream  of  the  flame  region  and  the  modified  impedances  ZQ  and  ZL  are 
given  by 


2  =  Z  +  M  and  Z.  =  Z.  +  M, 

OOO  L  L  L 


The  pressure  spectrum  Spp  is  obtained  by.  multiplying  Eq.(13)  by  its 
complex  conjugate  and  taking  an  ensemble  average.  The  result  is  (for  x  =  0) 


S»P(X*W‘!S  V*-°> 


where  S  is  the  spectrum  of  the  unsteady  heat  release  from  the  combustion 
qq 

region. 

The  natural  frequencies  of  the  duct  may  be  obtained  by  minimizing  the 
magnitude  of  the  denominator  of  the  Green's  function.  For  example,  consider 
the  case  of  a  rigid  walled  duct  with  no  convection  losses,  which  is 

A  A 

acoustically  rigid  at  x=0  (i.e.,  ZQ  =  «)  and  open  at  the  end  x  =  L  (i.e.,  ZL 
=  0).  In  r.i'is  case,  the  natural  frequencies  are  given  by  the  relation 


cos  (or  +  0)  +  [yjf-  -  1][  cos  a  cos  p  ]  -  0 


EXPERIMENTAL  EFFORTS 


The  experimental  set  up  used  is  shown  in  Fig.  2.  It  consists  of  a 

2 

rectangular  duct  three  meters  long.  It  has  a  7.5  x  5  cm  cross  section  and 
consists  of  an  inlet,  combustor  and  exhaust  sections.  The  inlet  contains  an 
injector  whose  face  is  made  of  sintered  stainless  steel.  A  mixture  of  propane 
and  air  is  introduced  into  the  set  up  via  the  injector. A  v-shaped  flame  is 
stabilized  in  the  combustor  section  on  a  0.8  mm  diameter  nichrome  wire.  The 
wire  is  attached  to  the  combustor  windows  at  half  combustor  height  and  it  is 
heated  electrically  to  improve  its  flame  holding  characteristics.  The  windows 
are  made  of  quartz  and  allow  optical  access  of  the  flame  zone. 

The  flow  approaching  the  stabilizing  wire  is  parallel  and  uniform  and  a 
fine  wire  mesh  grid  located  8  cm  upstream  of  the  stabilizing  wire  acts  as  a 
flame  arrestor  in  case  of  a  flashback.  The  cold  flow  Reynolds  number  is  kept 
below  10,000  in  the  interest  of  maintaining  a  disturbance  free  flow.  Typical 
cold  flow  velocities  considered  are  in  the  range  of  1-2  m/sec.  The  combustor 
walls  are  water  cooled  enabling  wall  mounted  pressure  transducers  to  be  used 
to  monitor  the  acoustic  pressure  field.  In  addition,  temperature  measurements 
are  carried  out  by  means  of  thermocouple  junctions. 

Capabilities  for  measuring  spontaneous  CH  species  radiation  from  the 

flame  have  also  been  developed.  The  concentrations  of  these  species  are  a 
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measure  ot  the  reaction  rate  and  ,  thus,  the  heat  release  rate  and  are 
useful  in  describing  the  unsteady  combustion  field.  Thus,  the  autospectrum  of 
the  unsteady  heat  release  rate  is  proportional  to  that  of  the  radiation  and 
the  following  relation  is  assumed  to  hold 


where  B  is  a  positive,  real  constant. 


The  radiation  signal  was  measured  from  the  entire  flame  to  obtain  the 
net  unsteady  heat  release  spectrum  from  the  combustion  zone.  The  autospectrum 
of  a  typical  measurement  is  shown  in  Fig.  3.  According  to  Eqs.  (15)  and  (17) 
this  spectrum  may  be  regarded  as  the  source  for  the  corresponding  pressure 
fluctuations.  Note  that  the  spectrum  is  dominated  by  low  frequency 
components.  Detailed  investigations  of  the  flame  region  are  presented  in  Ref. 
8  where  it  is  shown  that  this  pattern  occurs  due  to  the  dominating  influence 
of  large  scale  structures  in  the  flame  region. 

Measurements  of  the  radiation  signal  from  the  flame  were  carried  out 
under  several  different  conditions.  The  fuel  fraction  was  changed  between 
approximately  2.1  %  to  2.4  %.  Measurements  at  higher  fuel  fractions  could  not 
be  made  due  to  the  appearance  of  combustion  instability  in  the  system  which 
resulted  in  flame  flash  back  (see  Ref.  8).  Measured  flame  temperatures  were 
in  the  1250  -  1400°K  range.  Experiments  were  also  run  in  a  shortened  version 
of  the  set  up  shown  in  Fig.  2  by  removing  the  exhaust  section.  This  had  the 
effect  of  increasing  the  natural  frequencies  of  the  set  up.  In  addition,  some 
experiments  were  performed  in  the  presence  of  an  additional  flame  holder  in 
the  combustor  section.  This  served  to  modify  the  flow  field  in  the  flame 
region  and  addressed  the  question  whether  the  flow  field  influences  the  sound 
generation  process  in  a  direct  fashion  (i.e,  is  the  sound  field  changed 
appreciably  by  the  flow  field  even  without  changes  in  the  unsteady  heat 
release  rate).  It  will  be  discussed  shortly. 

To  carry  out  the  theoretical  calculations  of  the  pressure  spectra,  the 

A  * 

values  of  ZQ,  ZL,  the  wall  loss  factor  x  and  the  constant  B  ( Eq. ( 17 ) )  must  be 
known.  Impedance  measurements  in  cold  flow  (i.e.,  flow  without  combustion) 
tests  using  an  impedance  tube  technique  and  an  available  computer  program 


n 


indicated  that  the  set  up  behaved  to  a  good  approximation  as  an  acoustically 

*  * 

closed-open  duct.  Thus,  ZL  and  1/ZQ  were  both  set  to  zero.  The  wall  loss 
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factor  x  was  modelled,  following  Pierce,  as 

X  *  (18) 

where  Bj  is  a  constant,  indeper  t  of  the  frequency  and  depends  only  upon 
the  characteristics  of  the  duct  walls 

To  estimate  B1  and  B,  a  trial  and  error  technique  was  used.  A  sample 
pressure  and  corresponding  radiation  spectrum  were  chosen.  Using  Eqs.  (15) 
and  (17),  theoretical  pressure  spectra  corresponding  to  different  values  of 
Bj  and  B  were  calculated.  The  values  of  Bj  and  B  yielding  the  best  fit  with 
the  experimental  pressure  spectrum  were  determined.  Once  determined,  and  B 
are  known  for  all  the  tests  run  on  the  experimental  set  up  since  they  depend 
only  upon  the  duct  and  radiation  data  acquisition  system  characteristics 
(which  remained  same  in  all  the  tests). 

The  comparison  between  the  theoretical  and  experimental  pressure  spectra 
corresponding  to  the  radiation  spectrum  in  Fig.  (3)  is  shown  in  Fig. (4). 
Using  the  same  values  of  B  and  Bj  (recall  the  wall  loss  factor  ^s  defined  for 
unit  length  and  remains  the  same),  a  comparison  between  the  theoretical  and 
experimental  pressure  spectra  for  the  shortened  duct  case  is  shown  in  Fig.  5. 
The  unsteady  radiation  spectrum  measured  for  this  case  is  shown  in  Fig.  6. 
The  agreement  in  both  these  cases  and  for  all  the  other  cases  tried  (but  not 
shown  here  for  brevity)  is  very  good  and  serves  to  verify  the  theory.  Note 
that  the  comparison  is  made  only  in  the  low  frequency  range  (up  to 
approximately  250  Hz).  In  this  frequency  range,  the  flame  zone  length  (of  the 
order  of  20  cm)  is  an  order  of  magnitude  smaller  than  the  relevant  acoustic 
wavelengths. 


Some  comments  are  in  order  regarding  the  connection  between  the 

frequency  content  of  the  source  (  the  unsteady  heat  release  rate)  and  the 

pressure  fluctuations.  The  unsteady  heat  release  is  broadband  in  nature  (see 

Figures  3  and  6)  with  low  frequency  components  dominating  (frequencies  less 

than  approximately  100  Hz  in  the  considered  experiments).  As  shown  in  Ref.  8, 

this  behavior  is  due  to  the  presence  of  large  scale  structures  in  the  flame 

region.  For  an  unenclosed  flame,  this  frequency  content  is  reflected  in  the 
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resulting  pressure  fluctuations  which  also  are  broadband  in  nature  .  For 

enclosed  flames,  however,  such  as  those  considered  herein,  the  pressure 
spectrum  is  modified  by  the  presence  of  the  duct  natural  frequencies. 

Theoretically,  this  may  be  explained  in  terms  of  the  Green's  function  G  and 

the  requirement  that  the  frequency  spectrum  of  the  time  rate  of  variation  of 
the  heat  release  rate  (the  source  for  the  pressure  fluctuations,  see  Eq.  (7)) 
be  bounded.  The  latter  requirement  implies  that  (  varies 

approximately  as  w~n  where  n  is  a  positive  number.  This  implies  that  w  Sqq 
also  varies  as  uTn.  Thus,  the  frequency  dependence  of  the  pressure  spectrum 
may  be  given  by  (see  Eq.  (15)) 

c  ~n  |  r  t  2 

Spp  w  >  G  ■ 

As  noted  earlier,  the  Green's  function  G  is  maximized  at  the  natural 

frequencies  of  the  duct.  However,  the  presence  of  w  n  ensures  that  the 
pressure  fluctuations  also  exhibit  an  upper  cutoff  and  it  will  be  only  the 
first  few  natural  frequencies  that  will  be  significant  in  the  spectrum. 

The  theory  constructed  herein  is  valid  for  any  flame  shape  in  the 
combustion  region,  as  long  as  this  region  is  compact  compared  to  the  relevant 
acoustic  wavelength.  The  flame  shape  is  largely  dependent  upon  the  nature  of 
the  flow  f;nld  in  the  combustion  region.  Jn  the  theory,  the  flow  effects  are 


not  included  directly  but  are  included  implicitly  in  the  unsteady  heat 
release  (which  is  controlled  by  the  flow  to  a  large  extent)  which  is  presumed 
known.  Thus,  it  is  assumed  that  the  sound  pressure  fluctuations  are  mainly 
affected  by  changes  in  the  unsteady  heat  release  and  are  not  directly 
modified  by  the  flow  field  (i.e.,  the  effect  of  the  flow  field  is  indirect 
via  the  unsteady  heat  release).  To  check  this  assumption  ,  the  flame  shape 
was  changed  by  introducing  a  second  flame  holder  into  the  combustion  region 
as  shown  in  Fig.  7.  This  second  flame  holder  was  a  cylindrical  rod,  5  mm  in 
diameter  and  its  introduction  resulted  in  the  formation  of  a  double  v-shaped 
flame.  The  position  of  this  flame  holder  could  be  adjusted  with  respect  to 
the  original  flame  holder  (i.e.,  the  wire).  Thus,  several  different 
combustion  zones  could  be  investigated. 

As  an  example  consider  a  set  of  experiments  in  which  the  vertical 
displacement  Y  between  the  flame  holders  was  held  fixed  while  the  axial 
displacement  X  was  varied.  It  was  found  that  depending  upon  the  distance,  X, 
the  pressure  levels  measured  at  x  =  0  varied  between  100  to  140  db.  However, 
in  all  cases,  the  changes  in  the  pressure  levels  were  linearly  related  to  the 
changes  in  the  radiation  levels.  This  is  shown  in  Fig.  8,  which  plots  the 
obtained  pressure  levels  as  a  function  of  the  corresponding  radiation  levels 
at  the  first  natural  frequency.  The  predicted  pressure  spectra  were  also 
found  to  be  in  good  agreement  with  the  experimentally  measured  spectra  in  all 
cases.  These  experiments,  therefore,  verify  the  assumption  inherent  in  the 
theoretical  model,  that  changes  in  the  flow  field  affect  the  generated  sound 
mainly  by  changing  the  associated  combustion  rates  and  not  in  a  direct 


manner. 


CONCLUSION 


A  theoretical  model  capable  of  predicting  the  sound  generated  by 
confined  flames  has  been  developed  and  verified  by  comparison  with 
experimental  data.  The  relationship  between  the  frequency  contents  of  the 
unsteady  flame  heat  release  and  the  generated  pressure  spectra  has  been  shown 
to  depend  upon  the  natural  frequencies  of  the  duct  enclosing  the  flame. 
Finally,  it  has  been  shown  experimentally  that  for  low  velocity  situations 
the  flow  field  influences  the  sound  field  only  indirectly  through  its 
influence  upon  the  combustion  process. 
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